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ABSTRACT 

We describe a moderate-resolution FUSE survey of H2 along 70 sight lines to the Small and Large Magellanic 
Clouds, using hot stars as background sources. FUSE spectra of 67% of observed Magellanic Cloud sources 
(52% of LMC and 92% of SMC) exhibit absorption lines from the H2 Lyman and Werner bands between 912 and 
1120 A. Our survey is sensitive to N(H2) > 10'"^ cm"^; the highest column densities are log N(H2) = 19.9 in the 
LMC and 20.6 in the SMC. We find reduced H2 abundances in the Magellanic Clouds relative to the MiUcy Way, 
with average molecular fractions (fm) = 0.010!jj{jg^ for the SMC and (fm) = 0.012![J;[|JJf for the LMC, compared 
with (fm) = 0.095 for the Galactic disk over a similar range of reddening. The dominant uncertainty in this 
measurement results from the systematic differences between 21 cm radio emission and Lya in pencil-beam sight 
lines as measures of N(H I). These results imply that the diffuse H2 masses of the LMC and SMC are 8 x 10^ Mq 
and 2 x 10'' Mq, respectively, 2% and 0.5% of the H I masses derived from 21 cm emission measurements. The 
LMC and SMC abundance patterns can be reproduced in ensembles of model clouds with a reduced H2 formation 
rate coefficient, R ^3x 10"'^ cm-* s"', and incident radiation fields ranging from 10 - 100 times the Galactic mean 
value. We find that these high-radiation, low-formation-rate models can also explain the enhanced N(4)/N(2) and 
N(5)/N(3) rotational excitation ratios in the Clouds. We use H2 column densities in low rotational states (7 = 
and 1) to derive kinetic and/or rotational temperatures of diffuse interstellar gas, and find that the distribution of 
rotational temperatures is similar to Galactic gas, with (Tqi) = 82±21 K for clouds with N(H2) > lO""'^ cm"^. 
There is only a weak correlation between detected H2 and far-infrared fluxes as determined by IRAS, perhaps 
due to differences in the survey techniques. We find that the surface density of H2 probed by our pencil-beam 
sight lines is far lower than that predicted from the surface brightness of dust in IRAS maps. We discuss the 
implications of this work for theories of star formation in low-metallicity environments. 



L INTRODUCTION 

Molecular hydrogen (H2) must play a central role in our un- 
derstanding of interstellar chemistry, but little is known about 
the distribution of diffuse H2 in the interstellar medium (ISM) 
of the Galaxy. Major uncertainties remain about its formation, 
destruction, and recycling into dense clouds and stars (see the 
review by Shull & Beckwith 1982). Since the first detection of 
interstellar H2 (Carruthers 1970) numerous studies (Coperni- 
cus, Savage et al. 1977, hereafter S77; ORFEUS, Richter 2000; 
FUSE, Shull et al. 2000) using ultraviolet absorption measure- 
ments of H2 have created a general view of H2 formation, de- 
struction, and excitation in the diffuse ISM. However, we lack 
information about how the abundance of H2 and its physical 
parameters depend on the environmental conditions of interstel- 
lar gas, such as the metallicity, dust content, and UV radiation 
field. The H2 molecule has not been studied comprehensively 
in interstellar environments outside the Galactic disk. Varying 
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metallicity and UV radiation may modify the molecular abun- 
dance, thereby affecting the interstellar chemistry that triggers 
the star formation process. In particular, the formation rate of 
H2 in the diffuse ISM may depend on the composition and phys- 
ical state of the gas and grains in unknown fashion. 

Most of our knowledge of H2 in the diffuse Galactic ISM 
comes from ultraviolet absorption measurements of the Lyman 
and Werner rotational-vibrational bands in the far ultraviolet 
(912 - 1120 A). The Copernicus satellite mapped out the dis- 
tribution and properties of H2 along sight lines to stars confined 
to the local Galactic disk (Spitzer, Cochran, & Hirshfeld 1974; 
S77), and thus to a narrow range of gas properties. Despite 
this limitation, these studies uncovered important correlations 
between H2 abundance and the environmental conditions (dust 
and gas). The basic volume limitation on studies of H2 per- 
sisted until the advent of new instruments. The more sensi- 
tive but short-lived instruments, HUT (Gunderson, Clayton, & 
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Green 1998) and ORFEUS (Dixon, Hurwitz, & Bowyer 1998; 
Ryu et al. 2000; Richter 2000), have extended the range of H2 
studies to selected distant sources in the Galaxy and the Magel- 
lanic Clouds, and the IMAPS spectrograph (Jenkins & Peimbert 
1997; Jenkins et al. 2000) obtained high-resolution H2 spectra 
of bright targets in diverse environments. 

Molecular hydrogen is thought to form on interstellar dust 
grains when atoms are adsorbed on the grain surface, chemi- 
cally bond there, and then are ejected from the grain surface 
(Hollenbach, Werner, & Salpeter 1971). Theoretical calcula- 
tions, accounting for the probabilities of atom adsorption, and 
molecule formation and ejection, place the volume formation 
rate on grain surfaces at nanmR, where /? « 3 x 10"^' cm^ s"' 
in typical interstellar conditions, and where it is assumed that 
the total hydrogen density, is proportional to the number 
density of grains (Hollenbach & McKee 1979). Observations 
by Copernicus support this idea, with an inferred formation 
rate coefficient /? = (1 -3) x 10"" cm^ s"' (Jura 1974). Cor- 
relations of H2 and dust can provide a key test of this theory, 
but the narrow range of dust and gas properties in the Galactic 
disk accessible to Copernicus prohibited an exploration of the 
dependence of molecule formation on dust properties in a vari- 
ety of environments. The lower- me tallicity environments of the 
SMC and LMC (Welty et al. 1997; Welty et al. 1999) allow us 
to probe H2 formation and destruction in physical and chemical 
environments different from the Galaxy. In particular, the 4- 
to 17- times lower dust-to-gas ratios in the Magellanic Clouds 
(Koornneef 1982; Fitzpatrick 1985) imply a smaller grain sur- 
face area per hydrogen atom and a correspondingly lower effi- 
ciency of H2 formation. 

The Magellanic Clouds are therefore a valuable nearby lab- 
oratory for studying H2 in more distant, low-metallicity star 
forming regions and QSO absorption-line systems with high H 1 
column density. Molecular hydrogen has been detected in three 
damped Lya systems with N(H2)/N(H I) in the range 10"^ to 
IQ-'* (Ge & Bechtold 1999; Ge, Bechtold, & Kulkarni 2001; 
Petitjean, Srianand, & Ledoux 2000; Black, Chaffee, & Foltz 
1987), and limits have been placed on molecular abundance in 
several other systems. Some of these systems are known to 
have low metallicity and/or high UV radiation intensity relative 
to the Galaxy. Detailed studies of H2 in the nearby, spatially re- 
solved Magellanic Clouds can serve as an important benchmark 
for comparison with these high-redshift systems. 

Here, we describe the first survey results on H2 in the ISM of 
the Magellanic Clouds obtained with the Far Ultraviolet Spec- 
troscopic Explorer (FUSE) satellite. Initial FUSE results on H2 
in the Milky Way (ShuU et al. 2000) suggest that diffuse H2 is 
ubiquitous in the Galaxy. A large fraction of FUSE sight lines 
through the Galactic disk and halo exhibit absorption from the 
Lyman and Werner rotational-vibrational bands, showing that 
the Copernicus results on H2 in the local regions of the disk 
extend, in principle, to more distant regions of the Galaxy. The 
sensitivity of the FUSE spectrograph allows the use of distant 
hot stars and extragalactic objects as background sources, open- 
ing up more distant regions of the Galaxy and beyond to study 
of the H2 molecule. Our observations survey the H2 abundance 
and properties in the low-metalUcity regimes of the LMC and 
SMC, allowing us to model the formation and destruction rates 
of H2, and the density and radiation field in the interstellar gas. 

In § 2 we describe the observations and the data reduction 
and analysis schemes. In § 3 we survey the results on column 
densities, excitation, and other observed properties of the de- 
tected H2, discuss correlations with other gas and dust proper- 



ties along the sight lines, and compare these results to numerical 
cloud models. Section 4 discusses these results and summarizes 
the general conclusions of this survey. In this paper we assume 
a distance of 60 kpc to the SMC and 50 kpc to the LMC. 

2. OBSERVATIONS AND ANALYSIS 
2.1. FUSE Observations 

The observations reported here comprise all the available 
LMC and SMC targets from FUSE Guaranteed Time Obser- 
vations during Cycle 1 (up to October 2000). The FUSE mis- 
sion and its instrumental capabilities are described by Moos et 
al. (2000) and Sahnow et al. (2000). The target stars were se- 
lected for campaigns to study O, B, and Wolf-Rayet stars (Pro- 
gram ID PI 17) and to examine hot gas in the Milky Way (MW) 
and Magellanic Clouds (P103). Seven of the SMC stars was 
chosen specifically for H2 studies (PI 15). These observations 
present the first opportunity to study H2 throughout the Mag- 
ellanic Clouds in a comprehensive fashion. We include in the 
sample three SMC stars and one LMC star previously analyzed 
by Shull et al. (2000), the LMC star Sk -67 05, previously ana- 
lyzed by Friedman et al. (2000) and the SMC star Sk 108, pre- 
viously analyzed by Mallouris et al. (2001). An atlas of Mag- 
ellanic Clouds stars observed with FUSE has been compiled by 
Danforth et al. (2001). We list the target stars, dataset names, 
and observational parameters in Tables 1 and 2. A summary of 
results appears in Tables 3 and 4, where we list detected colunm 
densities N(H2), or upper limits, for the program stars. Tables 
5 and 6 list individual rotational level column densities N(J) for 
the LMC and SMC, respectively, as derived from the curve-of- 
growth and/or line-profile fits (§ 2.2). 

All observations were obtained in time-tag (TTAG) mode 
through the 30" x 30" (LWRS) apertures. Most observations 
were broken into multiple exposures taken over consecutive 
orbital viewing periods (see Tables 1 and 2). The photon 
lists from these exposures were concatenated before being pro- 
cessed by the current version of the calibration pipeline soft- 
ware (CALFUSE vl.8.7). This program computes the shifts in 
the detected position of a photon required to correct for: (a) 
the motion of the satellite; (b) nodding motions of the diffrac- 
tion gratings, which are induced by thermal variations on an 
orbital time scale; (c) small, thermally-induced drifts in the 
read-out electronics of the detector; and (d) fixed geometric dis- 
tortions in the detector. AppUcation of these shifts produces a 
two-dimensional, distortion-corrected image of a detector seg- 
ment, from which a small, uniform background is subtracted. 
We extract one-dimensional spectra for the LiF and SiC chan- 
nels recorded on each of the four detector segments. After cor- 
recting for detector dead time (which is small for count rates 
typical of targets in the Magellanic Clouds), we apply the most 
recent wavelength and effective-area calibrations to convert the 
detector count rate in pixel space to flux units as a function 
of heliocentric wavelength. Throughout these processing steps, 
we propagate 1 a uncertainties along with the data. No further 
modifications of the CALFUSE-processed data files are neces- 
sary for the H2 analysis. 

The FUSE detectors have very low dark count rates, and the 
overall level of light scattered to the detectors is low compared 
with the fluxes of our targets. Our targets range in specific flux 
from Fx = 5 X 10"^^ to 2 x 10"'^ erg cm"^ s"' A"'. The typical 
background count rate is the sum of dark counts and scattered 
light and corresponds to < lO"'"* erg cm"-^ s"^ A"'. The CAL- 
FUSE software performs a background subtraction, making the 
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Fig. 1 . — The product for the Lyman 0-0 to 9-0 bands of H2. Oscihator strengths and wavelengths are from Abgrall et al. (1993a,b). The best constraints on the H2 curve of growth 
are provided by the relatively weak 0-0 and 1-0 bands. The right axis shows the 4 cr column density sensitivity for the lines displayed, assuming that the limiting equivalent width of all 
lines is 30 mA (see Equation 1 and its discussion). 



contribution of the background to the total error in continuum 
placement and equivalent widths negligible. 

Because FUSE contains no internal wavelength-calibration 
source, all data are calibrated using a wavelength solution de- 
rived from in-orbit observations of sources with well-studied in- 
terstellar components (Sahnow et al. 2000). This process leads 
to relative wavelength errors of 10 - 20 km s"' across the band, 
in addition to a wavelength zero point that varies between ob- 
servations. Thus, we cannot rely completely on radial velocity 
measurements of interstellar absorption. In the analysis below, 
we derive relative velocities from H2 lines and fix the zero point 
in each sight line by assuming that the Milky Way absorption 
lies at Vlsr = km s"' . 

The resolution of the FUSE spectrograph across the band was 
A/AA ~ 10,000-20,000 for these observations. Because many 
of these observations were acquired during testing and calibra- 
tion phases of the FUSE mission, and because others were ob- 
tained late in 2000 after the instrument and its calibration had 
settled, the instrumental resolution and other parameters vary 
from target to target. However, since the results rely on mea- 
sured equivalent widths and damping-profile fitting, our conclu- 
sions are not sensitive to the resolution and spectrophotometric 
calibration. 



2.2. FUSE Data Analysis 

In our analysis of H2, we focus on the Lyman and Werner 
electronic transitions, with some 400 vibrational-rotational 
lines arising from the 7 = 0-7 rotational states of the ground 
vibrational and electronic states of H2. We search for 7 < 7, 
but lines above 7 = 4 are difficult to detect in diffuse interstellar 
clouds, given the typical 4 cr limiting equivalent width of 30 - 
40 mA and the corresponding column density limit of N(H2) ~ 
10'"^ cm"^. In this analysis we use oscillator strengths, wave- 
lengths, and damping constants from Abgrall et al. (1993a,b). 
The Hne strengths /A^ for the R(0), R(l) and P(l) Hnes of the 
0-0 to 9-0 Lyman bands appear in Figure Starting at 0-0, the 
line strengths increase by roughly a factor of ten to 7-0 and then 
decline slowly as the upper vibrational level increases. 

Figure ^ shows portions of FUSE spectra of four program 
stars, with a range of column densities N(H2). We show the 
Lyman (4-0) vibrational band and label Milky Way and Mag- 
ellanic lines. Absorption from gas in the Magellanic Clouds is 
easily distinguished by a large velocity separation in the range 
V = 200 - 300 km s"' for the LMC and 100 - 170 km s"' for 
the SMC. This paper is concerned with the Magellanic Cloud 
absorption only. Results on the detected Galactic gas in these 
sight lines will be presented later, together with other FUSE 
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Fig. 2. — Examples of FUSE spectra of SMC and LMC hot stars showing varying column densities of Ht. First panel: Sk -67 1 1 1 shows no detectable Ht at LMC velocity, with N(H2) 
< 5.0 X 10 cm"". The Galactic H2 lines are marked here for reference in the lower panels. Second panel: Sk -65 22, with N{H2) =1.6 X 10 cm"". Third panel: Sk -69 246, with N(H2) 
= 6.3x10 cm"". Fourth panel: The SMC star Sk 82, with N(H2) = 7.8 X 10 cm" . The wavelength range is chosen to encompass the Lyman 4-0 band of H2. The Lyman 5-0 P(4) line 
appears at left in the lower two panels. 



targets at high Galactic latitude. 

Figure || shows the spectra of three selected stars, plotted 
in velocity space to show the relative positions and strengths 
of the H2 lines of interest and the common blending between 
Galactic and Magellanic components. Figure ^ shows the en- 
tire FUSE spectrum of one LMC star, Sk -67 166. This figure 
illustrates the richness and complexity of FUSE data on H2. We 
have identified the Lyman (red) and Werner (green) bands and 
the MW and LMC components. The LMC H2 component has 
N(H2) = 5.5 X 10'^ cm"^ and the MW component has N(H2) = 
5.0 X lO'^cm-2. 

We employ a complement of techniques when analyzing H2 
absorption lines. The final products of the analysis are mea- 
surements of the column densities in the individual rotational 
levels, N(J), from which we can infer the gas density, radiation 
field, and formation and destruction rates of H2. For absorbers 



with N(H2) <^ 10 cm , we measure equivalent widths, W\, 
of all H2 lines and produce a curve of growth (COG) to infer 
a Doppler /7-parameter and N(7). For high-column density ab- 
sorbers with damping wings in the J = and 1 lines, we fit the 
line profiles to derive N(0) and N(l), and we use a curve-of- 
growth technique to derive column densities for J >2. How- 
ever, there is no well-defined column density limit above which 
profile-fitting must be performed. Instead, it is used in any case 
for which no lines of J = and 1 can be fitted with Gaussian 
line profiles, due to the presence of damping wings. 

We are often forced to neglect H2 lines in spectral regions 
near strong interstellar absorption or bright geocoronal emis- 
sion. The H2 Lyman (6-0) band is always lost in the strong 
damping wings of the interstellar Ly/3 line (1025.7 A), and the 
Lyman (5-0) band lies among the resonance absorption lines 
C II A1036.3 and C II* A1037.0. We neglect these bands ex- 
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cept when their J >2 lines appear apart from the intervening 
absorption. 

The key aspect of our analysis software is the rapid and con- 
sistent measurement of as many individual H2 absorption lines 
as possible in each sight line. Variations in the data quality, 
line-of-sight structure, and spectral type of the source prevent 
the use of a uniform set of H2 lines in the analysis. Thus, the 
process of measuring the numerous lines that enter the curve- 
of-growth fit cannot be automated completely. Our software re- 
quires the user to decide on a Une-by-line basis which H2 lines 
will be fitted. 

First, model spectra of varying N(J) and b are overlaid on 
the spectrum to obtain an estimate of the total column den- 
sity N(H2) and the radial velocities of the Galactic and Mag- 
ellanic components. At this stage, detections are discriminated 
from non-detections and analyzed separately. We describe our 
analysis scheme for detections here and return to discuss non- 
detections below. Our software scans for H2 lines from a list 
and queries the user at each line position. If the user decides the 
line is present and not severely compromised by blending, the 
line is fitted with a Gaussian profile, and its central wavelength 
(Ao), equivalent width (Wa), and full width at half maximum 
(FWHM) are stored in a table. The uncertainty in Wx is cal- 
culated, taking into account statistical errors in the data points 
and systematic uncertainty in the placement of the local contin- 
uum. A Une scan is performed separately for each FUSE de- 
tector segment to avoid introducing systematic errors that may 
be associated with combining the overlapping segments into 
one dataset. Segment-by-segment analysis is common practice 
among FUSE investigators and we adopt it here. After each 
scan, lines that are blended or that were excluded from the scan 
for some other reason are fitted separately and placed in the 
table. There is one table for each of the four FUSE detector 
segments (la,lb,2a,2b) and two channels (LiF and SiC), for a 
total of eight files. 

At the conclusion of scanning, these tables contain the mea- 
surements for all the suitable H2 lines in the spectrum, includ- 
ing multiple measurements for those lines that appear on more 
than one detector segment. There is a maximum of four mea- 
surements for each line. Prior to COG fitting, these tables are 
merged with the following scheme. The final W\ measurement 
for each line is the average of the individual measurements, 
weighted by the inverse squares of their individual uncertain- 
ties. The uncertainty in the final measurement is the maximum 
of these three quantities: (1) the quantity (X)^?)"^' where the 
(T,- are the individual uncertainties; (2) one-half the difference 
between the smallest and largest measurements (as a proxy for 
the standard deviation); or (3) 10% of the weighted mean itself. 
For most lines, (1) obtains the maximum value. However, for 
weak lines (2) is useful if duplicate measurements are widely 
separated and the formal standard deviation has little meaning. 
Finally, (3) is a conservative assumption that attempts to ac- 
count for unquantified or unknown systematic errors in a uni- 
form fashion. 

After the individual line tables are combined, the software 
produces an automated fit to a curve of growth with a single 
Doppler b parameter. This simple assumption is maintained 
throughout and is never clearly violated. Tests performed on 
high-quality data could not confirm the need for multi-valued h 
increasing with J, an effect seen in Copernicus data (Spitzer et 
al. 1974) and attributed to unresolved components with differ- 
ent rotational excitation. We use a downhill- simplex method to 



minimize the reduced statistic, with column densities N(J) 
and b as parameters. We typically obtain reduced in the 
range 0.5 - 2.0, with lower values for higher signal-to-noise 
data. 

The simplex -minimization technique is fully automated 
and well suited to deriving best-fit column densities and b, but 
it is inadequate for deriving realistic uncertainties in the mea- 
sured parameters. The usual scheme of using the covariance 
matrix at the best-fit point to derive confidence intervals on 
the parameters is not appropriate in the highly nonlinear, multi- 
dimensional parameter space of the curve of growth. Typically, 
the quantity d^^ /dpi varies substantially with parameter p,, 
such its value at the best-fit point does not accurately describe 
a midpoint between the critical values corresponding to the 
1 (J confidence interval. So, while we use a standard minimiza- 
tion technique for finding the best-fit column densities and b, 
we take a different approach to finding the maximum Ict varia- 
tions in the parameters (the error bars in Tables 5 and 6). 

In this multi-dimensional parameter space, the column densi- 
ties N(J) and the Doppler b parameter have somewhat different 
status, with the uncertainty in b contributing most of the uncer- 
tainty in N(J). We first calculate the critical value, x?rit' corre- 
sponding to the 1 a confidence interval. This quantity is given 
by the best-fit x^ plus an offset Ax^, which is distributed for 

parameters like x^ for N degrees of freedom (Bevington & 
Robinson 1992). Over a range of values containing the best 
fit, the software produces a curve-of-growth fit over a succes- 
sion of b values, in 0.1 km s~' intervals. At each fixed b, the 
best fit is obtained by allowing N(J) to vary freely. If this best 
fit has x^ < Xciit' then the column densities for different J are 
varied randomly to explore the parameter space in the immedi- 
ate region of the best fit with fixed b. Randomly chosen column 
densities that satisfy the x^ criterion are stored in a table for 
further analysis. This process is repeated over a range of fixed 
b values, building a table of permitted parameters that satisfy 
the x^ criterion. The quoted uncertainties in the published N(J) 
and b represent the maximum variation among all the parame- 
ter sets that have x^ < Xcrif '^^'^ '^ey represent the true 1 cr joint 
confidence interval on the parameters. 

In addition to its efficiency, this technique has the further 
benefit of providing realistic asymmetric error bars on the mea- 
surements. This technique is more efficient than a determin- 
istic multi-dimensional grid, which requires many more COG 
evaluations if it is "square" in the parameters. By contrast, our 
technique does not vary the column densities if the best fit for 
a given b is excluded. It requires approximately 10,000 COG 
evaluations for each sight line, far less than a deterministic grid 
with sufficient resolution to capture all the variations in the pa- 
rameter space. 

After COG fitting and exhaustive parameter space explo- 
ration, we have derived N(J), b, and uncertainties for each line 
of sight. However, for diffuse clouds with N(H2) ^ 10'** cm"^, 
it is better to derive N(0) and N(l) by fitting the wings of the 
damping profiles. For fines on the damping part of the COG, we 
combine the COG fitting technique described above with profile 
fitting. Model line profiles are produced with parameters N(J) 
and b for each H2 component in the spectrum, and convolved by 
a Gaussian instrumental line spread function. Typically, there 
are two components, corresponding to Galactic and Magellanic 
gas. Fitting is performed over the range of wavelength neces- 
sary to capture all the relevant lines and exclude difficult regions 
of the stellar continuum. The local continuum in the region of 
the fitted line profiles is described by a quadratic function. We 
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also add atomic and ionic lines where necessary to improve the 
fit. Although H2 lines with J > 2 are included in the fitting, their 
final column densities are instead calculated using the COG fit- 
ting technique. Finally, we fit the Lyman bands separately on 
each FUSE detector segment and combine the measurements in 
a weighted average for final results. The uncertainties quoted 
in Tables 3-6 are the standard deviations of the multiple mea- 
surements (a more thorough discussion of this fitting is given by 
Rachford et al. 2001). We find remarkable consistency among 
the different Lyman bands and the different detector segments, 
giving us confidence that 10% uncertainty can be achieved rou- 
tinely with this fitting method. 

In a few cases, we derive very small uncertainties on N(0) 
and N(l). The most extreme case is Sk -69 246 for which the 
uncertainties are less than 10%. To test whether the data quality 
supports such small uncertainties, we have performed a Monte 
Carlo simulation of the effects of noise on the fits. We generate 
a noiseless profile based on the observed column densities, and 
then add a series of random noise vectors whose characteristics 
match the observed noise in the object spectra. For Sk -69 246 
we generated 10 noisy synthetic spectra and performed profile 
fits. The standard deviations of the values of N(0) and N(l) 
were 0.01 dex, smaller than the reported uncertainties. This re- 
sult supports the claim that the segment-to-segment variations 
in the fitted values dominate the total uncertainty. 

In our line analysis, COG fitting, and profile fitting, we have 
paid careful attention to statistical and systematic uncertain- 
ties in deriving H2 column densities and Doppler parameters. 
The major source of uncertainty in the column densities derived 
from COG fitting is the often weak constraint on the Doppler 
b parameter Generally, lines from a given rotational level with 
a wide range of line strengths (/A, the oscillator strength times 
the rest wavelength in A) are necessary to adequately constrain 
b for a given sight line (see Figure [l|). In practice, the best con- 
straints are provided by the relatively weak Lyman 1 -0 and 0-0 
bands. If lines with both high and low /A are available, or for 
lines on the linear portion of the COG, the logarithmic errors in 
the column densities can approach 0.10 dex, and vary with the 
signal-to-noise ratio. If the available lines for a given J-level are 
spaced narrowly in /A, the error in this column density will be 
determined largely by the permissible range in b. In this case, 
solutions for the N(J) are coupled together by the b value; low- 
ering or raising b will move all the N(J) up or down together, 
respectively. Large uncertainties in b can also cause very asym- 
metric error bars on the column densities. If the best-fit column 
density for a given b lies near the best-fit N(J) for the extreme 
end of the b range, then the N(J) may take on a large uncertainty 
in the other direction, corresponding to the other extreme of b. 
All these systematic variations make the ratios in the column 
densities (seen below in § 3) less uncertain than the individual 
column densities themselves. 

A further systematic uncertainty is due to the unknown com- 
ponent structure of the H2 clouds. At the 15 - 30 km s"' res- 
olution of FUSE, we generally cannot resolve the fine-scale 
structure of complicated interstellar components observed in 
atomic lines with higher resolution (e.g., by Welty et al. 1997, 
1999). Our analysis assumes that all the H2 resides in one 
component, and that absorption from all rotational levels can 
be described by a single Doppler b. Snow et al. (2000) and 
Rachford et al. (2001) used the observed Na I component struc- 
ture along Galactic translucent-cloud sight lines to construct a 
multicomponent curve of growth. Here, the stronger absorp- 



tion makes it more likely that H2 is contained in more than one 
component. They found that, while the measured column den- 
sities N(J) were sensitive to the component structure, the dif- 
ferences in N(J) for lines on the Doppler (flat) portion of the 
COG were small compared to the uncertainties associated with 
the flat COG. Thus, we assume single-component structure and 
leave more detailed analysis for future examination of individ- 
ual sight lines from our survey. 

The sight lines in which H2 is not detected receive a separate 
analysis. If there is no visible evidence of the strong Lyman 
7-0 R(0) and R(l) lines, we place upper limits on the equiv- 
alent widths of these lines and convert these limits to column 
density limits assuming a linear curve of growth. We use the 
following expression for the (4 cr) limiting equivalent width of 
an unresolved line at wavelength Aq: 



(A/AA)(5/A^)' 

where A/ A A is the spectral resolution, and S/N refers to the 
signal-to-noise ratio per resolution element in the 1-2 A region 
surrounding the expected line position. By contrast, in Tables 1 
and 2 we list the S/N ratio per pixel to provide a means of com- 
paring datasets in a resolution-independent fashion. Because 
the resolution of FUSE varies across the band and even between 
observations, we conservatively assume A/AA = 10,000 for all 
upper limits. The N(0) and N(l) limits in Tables 5 and 6 are 
imposed with the strongest J = and J = 1 Lyman lines, 7-0 
R(0) and 7-0 R(l). For a linear COG, these lines have equiva- 
lent widths Wx = (26.8 mA) (N(0) / 10'"* cm'^) and Wx = (18.7 
mA) (N(l) / lO''* cm"^), respectively. The entries in Tables 3 
and 4 show a limit on the sum of the two individual lines. Most 
of the upper limits lie just below the lowest detections of H2. 

In summary, we describe the results and uncertainties that 
appear in Tables 5 and 6. We list the final column density 
and uncertainty results for aU our program stars. These col- 
umn densities were derived from the COG or profile-fitting 
routines described above. The uncertainties were calculated 
from either COG parameter- space searches or band-to-band 
and segment-to-segment variations in the profile fits, again as 
described above. The 4 a upper limits are derived according to 
Equation (1) from individual signal-to-noise ratios of the data 
in the region of of the Lyman 7-0 R(0) and R(l) lines. For these 
limits we assume an effective resolution R = 10,000. For low- 
column density sight lines where no constraint on the Doppler 
b parameter is provided, we assume that the lines fall on the 
linear portion of the curve of growth and label these cases ac- 
cordingly. These tables form the fundamental database for all 
subsequent analysis. 

2.3. Notes on Individual Sight Lines 

The LMC targets Sk -69 243 (R136) and MK 42 He at the 
center of the 30 Doradus H II region, and Sk -69 246 lies in an 
isolated field 4' to the north. The Sk -69 246 spectrum shows 
no evidence of unusual properties. The other two datasets, how- 
ever, show line profiles much broader than the expected FUSE 
line spread function in both the LMC and Galactic absorption 
lines due to crowding among the 70h- hot stars in the inner 30" 
of the cluster. Because the blurring mimics a low-resolution 
spectrograph, we can only be sure that there is H2 at LMC 
velocities in these sight lines, with N(H2) ~ 10'^ cm"^. The 
crowded fields and line broadening preclude a detailed analysis 
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of the column densities. We include these stars in the H2 detec- 
tion statistics (Tables 1 and 3) but exclude them from the LMC 
column density sample (Table 5). 

The sight lines to NGC 346 (3, 4 and 6) lie within a 50" 
circle on the sky in the SMC. As seen in Tables 4 and 6, their 
column densities and Doppler b parameters are statistically in- 
distinguishable, suggesting that these closely spaced sight lines 
probe the same interstellar cloud. In the abundance and excita- 
tion distributions and tests performed below, we use the mean 
column density in each level and the mean b to stand in for these 
three sight lines. 

Mallouris et al. (2001) reported tentative detections of J = 1 
and J = 3 lines in towards the SMC star Sk 108. Our survey im- 
poses formal 4 a limits on the equivalent widths of undetected 
H2 lines. According to this criterion, these lines are not signif- 
icant, and we place an upper limit on N(H2) for Sk 108 that is 
larger than the value they report. 

For closely paired sight lines we can use the observed H2 ex- 
citation to estimate the sizes of diffuse interstellar clouds. The 
sight lines to HD 5980 and AV 232 are separated by 58" on the 
sky, and lie ~ 17 pc apart at the distance of the SMC. Shull et al. 
(2000) reported different rotational excitation N(J) in these two 
sight lines, indicating that they may probe different molecular 
gas. The Galactic components in these sight lines show sim- 
ilar column densities and kinetic temperatures, and are likely 
to arise from the same cloud. However, the SMC H2 compo- 
nents on these two sight lines show distinct differences in rota- 
tional excitation. Measurement of H2 excitation and abundance 
in multiply-intersected absorbers will allow us to constrain the 
sizes and structure of diffuse interstellar clouds. 

2.4. Ancillary Data 

To fully exploit the FUSE data on H2, we require additional 
information about the dust and gas along the sight Unes. The 
most important of the supplemental data is the atomic hydro- 
gen column density along each sight Une. To obtain N(H I) for 
our entire sample, we use 21 cm radio emission and apply a 
correction to account for known systematic errors. 

We obtain H I column densities from the recent Parkes and 
Australia Telescope Compact Array maps (SMC, Stanimirovic 
et al. 1999; LMC, Staveley-Smith et al. 2001) and then apply 
a correction, discussed below. These maps were created with 
60" (LMC) and 98" (SMC) beams and were provided to us 
in fully caUbrated, machine-readable form (L. Staveley-Smith, 
personal communication). To obtain the N(H I) for the FUSE 
sight lines, we select the four beam positions that surround each 
star and perform a bilinear interpolation. Using N(H I) from 21 
cm emission may be subject to systematic uncertainties associ- 
ated with dilution in the 1' beam if the interstellar gas is highly 
clumped on smaller scales or if there is additional emission be- 
hind the stars. Available UV absorption measurements provide 
a valuable check on these measurements. 

Attempts to derive N(H I) from the Ly/J lines in the FUSE 
band generally fail, owing to the severe blending between the 
Galactic and SMC/LMC components, geocoronal O I emission, 
and the presence of the H2 absorption. The higher Lyman series 
is present in the FUSE band but is severely compromised by the 
dense forest of H2 and atomic lines below 1000 A. Data on HI 
Lya are available for roughly two-thirds of our sample, but we 
find that reliable determinations of N(H 1) from a combination 
of Lya and Ly/? profile fitting is impractical in all but the best 
cases. The largest uncertainties are in the placement of the stel- 
lar continuum near these lines and in the close velocity spacing 



of the Galactic and Magellanic components (100 - 300 km s~^). 
While determinations could be made for roughly 30% of the 
targets with high-resolution, high-S/N data (from HST/FOS, 
GHRS, and particularly STIS) and simple continuum structure, 
we could not achieve this goal for the entire sample. We per- 
formed simple fits to the Lya profiles to estimate N(H 1) for 
the one-third of the sample where the profile is not severely 
compromised. In these trials, we found that the constraints on 
N(H I) range from 0.3 to 1.0 times the 21 cm value, with a 
mean of 0.56 for both Clouds. We have not attempted to re- 
move contaminating absorption and complicated features in the 
stellar continuum, which are present even in these best cases. 
Thus, fitting the Lya profile is itself uncertain, and generally 
establishes an upper limit. The lower limit is more difficult to 
determine because the Magellanic component can be reduced 
and the Galactic component increased to compensate. How- 
ever, it appears that for 25 sight lines, the upper limit derived 
from Lya is, on average, half its value from 21 cm emission. 
This result suggests that we should adjust our H I column den- 
sities to account for this systematic error. 

We are reluctant to apply a uniform factor of 0.5 correction to 
the entire FUSE sample, because the subsample may have inter- 
nal biases. For example, the stars with HST data may be located 
preferentially on the front sides of the Clouds. As a compro- 
mise between uniformity and uncertainty, we reduce the 21 cm 
N(H I) by a factor 0.75. We then impose a large uncertainty on 
N(H I) to account for this systematic error. For the final col- 
umn density, we adopt N(H 1) = (0.75![!;|) x N(21 cm). This 
scheme captures the uncertainty in the column density and can 
be applied to the entire sample. Both extremes, at the high end 
from 2 1 cm data and at the low end from our HST estimates, are 
included in the error bars. The final uncertainty in the molecular 
fraction is dominated by the uncertainty in N(H 1). 

The color excess, E(B- V), is an important measure of the dust 
abundance towards the target stars. The color excesses were ob- 
tained using two methods: (1) the observed B-V color and an 
adopted intrinsic (B-V)o scale (Fitzgerald 1970 and Schmidt- 
Kaler 1982, with -0.31 or -0.32 adopted for early giants and su- 
pergiants), and (2) the observed UV spectrophotometry (lUE) 
plus B and V together with model atmospheres to fit the spec- 
tral energy distribution. In practice, we take the mean of the 
two methods to minimize any external uncertainties. Because 
intrinsic colors for Wolf-Rayet (WR) stars vary dramatically, 
model atmospheres together with UV and optical spectropho- 
tometry are necessary to derive E(B-V). We estimate the final 
uncertainty in these values to be < 0.02 mag. A further uncer- 
tainty is the contribution of foreground dust to the total sight 
line reddening. We derive the corrected color excess, E'(B-V), 
by subtracting from the measured color excess a foreground 
E(B-V) = 0.075 for the LMC and 0.037 for the SMC. These 
average values were obtained by the COBE/DIRBE project by 
averaging the dust emission in annuh surrounding the Clouds 
(Schlegel, Finkbeiner, & Davis 1998). In § 3 we describe a test 
of these averages with the observed gas-to-dust ratio. 

We use the Copernicus survey of Savage et al. (1977) to com- 
pare the similarities and differences between the FUSE sample 
and H2 in the disk of the Galaxy. The sample of 109 stars ob- 
tained by S77 (their Table 1) was culled for the 51 stars with 
measured N(H 1), either a measurement or limit on H2, and E(B- 
V) < 0.20. The 17 upper limits in this sample are included in 
certain tests depending on the context of the comparison. In ad- 
dition, we compiled column densities for y > 2 for 3 1 of these 
stars from Spitzer, Cochran, & Hirshfeld (1974). These column 
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densities are compared to the rotational excitation of the LMC 
and SMC gas in § 3. 

2.5. Sample Selection and Bias 

Most of the FUSE sample of LMC and SMC stars was not 
chosen specifically for the H2 survey. Most of the sample is 
drawn from three large FUSE team programs devoted to stel- 
lar and interstellar research. The stellar wind (Program PI 17) 
targets were chosen to occupy a range of spectral types. The al- 
lotted time for each observation was 4 ksec, so this is essentially 
a flux-limited sample. In practice, this criterion selects against 
stars with substantial extinction, which would be favored in a 
sample chosen to study H2. The stars in Program P103, the 
FUSE study of hot gas in the Magellanic Clouds, were cho- 
sen based on the availability of supporting observations and to 
probe the supershells and bubbles in the ISM of the LMC. In 
effect, this program selected against denser regions and high 
extinction in favor of regions where cavities of hot gas are ex- 
pected. In the SMC, seven stars were chosen specifically for 
H2 studies (Program PI 15). All others were drawn from the hot 
star or hot gas programs. However, from the point of view of 
this survey, the stars are located randomly across the LMC and 
SMC, but probably preferentially situated on the Milky Way 
side of the dense absorbing regions in the two galaxies. 

In addition to these selection biases, there is the additional 
problem that flux-limited or extinction-limited surveys favor 
targets on the foreground face of the LMC and SMC. This bias 
complicates the use of H I column densities derived from emis- 
sion measurements. The 21 cm emission beam samples all the 
gas, including some that lies behind the target star but does not 
appear in absorption along the FUSE pencil-beam sight lines. 
In this case, the measured gas-to-dust ratio, or N(H I) / E'(B- V), 
will be enhanced over the average ratio determined using H I 
Lyman series absorption to derive N(H I). Figure || compares 
the corrected color excesses for the program stars with their H I 
column densities from 21 cm emission, corrected as in § 2.4. 
The values of N(H I) / E'(B-V) range from lO^^- 10^^ cm'^ 
mag"' for the LMC stars, with a general trend towards decreas- 
ing values of the ratio with higher E'(B-V). There is a similar 
trend in the SMC, with N(H I) / E'(B-V) = lO^^ ^ - lQ2i-^ cmT^ 
mag"'. We plot the average value of this ratio for targets in the 
disk of the Milky Way as determined from lUE measurements 
by Shull & van Steenberg (1985), and the average values for the 
LMC (Koornneef 1982) and the SMC (Fitzpatrick 1985). Our 
target stars lie well above the Galactic average, but are consis- 
tent with the LMC and SMC values. Because E'(B-V) appears 
in both axes of this plot, the error bars on the points should run 
diagonally from the upper left to the lower right of each point. 
This fact probably explains the trend seen in the points, with 
lower E(B-V) lying above the average and higher E(B-V) be- 
low. We conclude that, while the location of the stars in the gas 
layers of the LMC and SMC may account for some of the scat- 
ter in the molecular fraction results, there is no large systematic 
increase in N(H I) and no corresponding systematic change in 
quantities derived from it. The observed trends suggest that our 
adopted correction to the 21 cm H I columns is reasonable. Had 
we adopted a larger correction to the 21 cm columns, the points 
that lie below the average gas-to-dust ratios would be far more 
discrepant. For further discussions of uncertainty in N(H I), see 
§2.4. 

3. SURVEY RESULTS AND ANALYSIS 



3.1. General Results 

A summary of results appears in Tables 3 and 4, where we 
list detected column densities N(H2), or upper limits, for the 
program stars. Tables 5 and 6 list individual rotational level 
column densities N(J) for the LMC and SMC, respectively, as 
derived from the curve-of-growth and/or line-profile fits (§ 2). 

The immediate, striking result of this survey, apparent from 
Tables 3-6, is the lower frequency of H2 detection in the LMC 
compared with the SMC. Shull et al. (2000) reported seven 
sight lines through the Galactic disk and halo, only one of 
which did not show detectable H2 (PKS 2155-304). The de- 
tection rate in the initial FUSE sample of extragalactic targets 
exceeds 90%. However, the total rate of detection in the LMC 
is 23 detections in 44 sight lines, for a 52% success rate. In the 
SMC, the detection rate is 24 of 26, or 92%. This difference 
appears to be related to the patchy structure of the LMC ISM, 
which contains many identified H I superbubbles blown out by 
OB associations and supernova remnants (Kim et al. 1999). By 
comparison, the ISM of the SMC is more quiescent and less 
patchy (see § 4-5 for further discussion). 

3.2. H2 Abundance 

The lower abundance of H2 in the low-metallicity gas of the 
LMC and SMC is a key result of this survey. The molecular 
hydrogen abundance is sensitive to the formation/destruction 
equilibrium of H2 in diffuse gas, where far-ultraviolet (FUV) 
radiation competes with formation on dust grains to determine 
the relative abundance of H2. The size of the sample allows 
us to apply statistical tests of H2 abundance in the LMC and 
SMC and compare it to H I, dust properties, and diagnostics 
of the UV radiation field. The reduced molecular fraction may 
indicate less efficient formation of H2 molecules on interstel- 
lar dust grains or enhanced H2 photo-dissociation by an intense 
FUV radiation field, or a combination of these two effects. In 
this section, we use simple analytic expressions and numerical 
models to assess the effects of varying environmental condi- 
tions on the H2 formation/destruction equilibrium. 

To assess the total abundance of H2, we define the quantity 
/h2 = 2N(H2) / [N(H I) + 2N(H2)], which expresses the fraction 
of hydrogen nuclei bound into H2. The molecular fractions for 
our sight lines appear in Figure ^ correlated with the corrected 
color excesses of the target stars. In Figure ^ we compare the 
S77 sample of Galactic disk stars with the LMC and SMC tar- 
gets. The key difference between the samples is the presence of 
several LMC and SMC sight fines with fm > lO""* at low E'(B- 
V). In the Clouds we find a wide range of molecular fractions at 
low reddening, E'(B-V) < 0.08. In the LMC there is a notice- 
able gap near log fm = -4, which does not appear in the SMC. 
The most striking pattern in these plots is the large number of 
sight lines to the LMC that show no detectable H2. S77 found 
low values of fm at low E(B-V) and a transition from low to 
high values { fm ~ 0.1) at a color excess E(B-V) = 0.08. This 
pattern is not readily discernible in the samples. 

Our results reverse prior indications that the average molec- 
ular fraction in the LMC is not reduced relative to the Galaxy. 
Using HUT, Gunderson et al. (1998) measured N(H2) ~ 10^° 
cm"2 for two LMC stars, Sk -66 19 and Sk -69 270, that do 
not appear in our sample. These sight lines have fm = 0.04, 
similar to the highest-/H2 points in the FUSE sample. They 
concluded that, in the low-metalHcity LMC, N(H2)/E(B-V) did 
not differ from the average Galactic value, despite large differ- 
ences in N(H I)/E(B-V). With the larger FUSE sample, we find 
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a substantial reduction in fm in the LMC. With ORFEUS data 
on three LMC stai's, Richter (2000) found that, although N(H2) 
relative to E'(B-V) in the Clouds might be similar to the Galac- 
tic values, the N(H2)/N(H I) ratio was lower than the Galactic 
average. The larger FUSE sample confirms this result. 

To quantitatively compare the distributions of fui for the 
Galaxy and Magellanic Clouds, we performed Kolmogorov- 
Smirnov (KS) tests on the samples. The two-sided KS test finds 
a 0.224 probability that the S77 and SMC samples are drawn 
from the same parent distribution, and a <0.001 probability for 
the S77 and LMC samples. Compared to each other, the SMC 
and LMC samples give a 0.121 probability of being drawn from 
the same distribution. This difference reflects the very different 
detection rates in the Clouds, but they suggest that the reduced 
fm in both the LMC and SMC is significant. 

For a quantitative comparison of the Galactic disk and LMC 
molecular fractions, we define the average column-density 
weighted molecular fraction (/) to be: 

(/) = J2 2N(H2)/ [N(H I) + 2N(H2)] . (2) 

In evaluating the average molecular fractions, we include the 
FUSE upper limits as measurements, so that (/) from FUSE is 
also an upper limit. For the LMC, we find (/) = 0.012 for stars 
with E'(B-V) = 0.00 - 0.20. For the SMC, we find (/) = 0.010 
for stars with E'(B-V) = 0.00 - 0.12. For the Galactic compar- 
ison sample, we find (/) = 0.039 for stars with E(B-V) = 0.00 
- 0.12, for comparison with the SMC. For Galactic stars with 
E(B-V) = 0.00 - 0.20, we find (/) = 0.095, for comparison with 
the LMC. If we exclude the upper limits in Tables 5 and 6 from 
the sample, we find (/) = 0.010 for the SMC and (/) = 0.021 
for the LMC, due to the large large number of non-detections 
there. 

Figure [s] shows the individual and cumulative molecular frac- 
tions correlated with the total H content, N(H) = N(H I) + 
2N(H2). The S77 sample shows a clear break near log N(H) 
= 20.7, above which all clouds have molecular fraction fui 
^0.1. This break is usually identified as the point at which the 
cloud achieves complete self-shielding from interstellar radia- 
tion (S77). The LMC and SMC sight lines show no such clear 
transition, even at the higher N(H) they probe. As we show be- 
low, the LMC and SMC show evidence of reduced formation 
rate of H2 on dust grains and of enhanced photo-dissociating 
radiation relative to the Galaxy. In this case, we expect to see 
a transition from small fm with large scatter to a tighter band 
of high /h2, but at a higher total H column density than in the 
Galactic disk. The FUSE samples require that such a transition 
occurs at logN(H) > 21.3 in the LMC and logN{H) > 22 in 
the SMC. Planned observations of several LMC and SMC stars 
with E'(B-V) > 0.20 will address this point in the future. 

The lower observed (/) is evidence for lower molecular 
abundance in the LMC and SMC. However, there is a poten- 
tial systematic effect that we must reemphasize. The S77 sur- 
vey used column densities, N(H I) and N(H2), determined from 
the same Copernicus spectra. This technique does not suf- 
fer the systematic uncertainty of comparing measurements of 
N(H I) from radio emission measurements with a large beam 
and N(H2) determined from pencil-beam sight lines. Thus, 
relative to the S77 measurements, we may underestimate the 
molecular fraction due to the systematic overestimate of the H I 
column determined from 21 cm emission if a substantial frac- 
tion of the H I seen in emission lies behind the target star or if 
clumping is significant below 1' scales. However, because the 



average dust-to-gas ratio for our targets scatter around the av- 
erage values found for the LMC and SMC, this problem is un- 
likely to exceed the generous error budget assigned to N(H I). 
(See the discussion of these issues in § 2). 

A simple model of H2 in an interstellar cloud (Jura 1975a,b) 
can be expressed as an equilibrium between the formation of H2 
with a rate coefficient R (in cm^ s~') and the photodestruction 
of H2 with rate D (in s"'): 

RnHiriHi = Dn(H2) ~ /d,ss ^ PiJ)n(H2,J), (3) 
J 

where P(J) is the photoabsorption rate from rotational level J 
and hh = n{H /) + 2n(H2). The coefficient /diss corresponds to 
the fraction of photoabsorptions that decay to the dissociating 
continuum and varies from 0. 10 to 0. 15 depending on the shape 
of the UV radiation spectrum. Molecules are destroyed by pho- 
toabsorption in the Lyman and Werner bands, followed by ra- 
diative decay to the vibrational continuum of the ground elec- 
tronic state. The formation of H2 is believed to occur on the 
surfaces of dust grains, but may also occur in the gas phase at a 
lower rate (Jenkins & Peimbert 1997). 

If we further assume a homogeneous cloud and replace the 
number densities «hi and n(H2) with column densities N(H I) 
and N(H2), respectively, then: 



where fm is the molecular fraction defined above. Thus, the ob- 
served molecular fraction is a measure of interaction between 
formation and destruction, and the molecular fraction can be 
suppressed both by reduced formation rate on dust grains and 
by enhanced photodissociating radiation. 

To assess the reduced molecular fraction, and, later, the ro- 
tational excitation of the sight lines, we use the results of an 
extensive grid of models produced by a new numerical code 
that solves the radiative transfer equations in the H2 lines and 
produces model column densities for a range of model parame- 
ters (Browning, Tumlinson, & Shull 2001, in preparation). This 
code has been developed with the FUSE surveys of H2 in mind. 
The model clouds are assumed to be one-dimensional and are 
illuminated on one side by a uniform radiation field of spec- 
ified mean intensity, /, in units photons cm"^ s"' Hz"'. (The 
code can illuminate the cloud on both sides, but the resulting 
changes in the model column densities are small.) generations 
of the code will permit two-sided clouds illuminated on both 
sides.) Because it follows the radiation field and H2 abundance 
in detail, the code makes no distinction between the optically 
thin and optically thick (self-shielded) regimes treated sepa- 
rately by Jura (1975a,b). The cloud density Wh, kinetic tem- 
perature Tkin, Doppler b parameter (to approximate thermal or 
turbulent broadening), cloud size d, and H2 formation rate coef- 
ficient on grains {R, in units cm"* s"' ) are the model parameters. 
The comparison grid calculated for this work consists of 3780 
models with / = (1.0, 4.0, 10, 20, 50, and 100) x 10"^ photons 
cm-2 s-i Hz-i, R = (0.1, 1.0, 3.0) xlO"'^ cm^ s"', = (5, 25, 
50, 100, 200, 400, 800) cm'^, d = (2, 4, 6, 8, 10) pc, and = 
(10, 30, 60, 90, 120, 150) K. For comparison with the LMC and 
SMC, we assume that "Galactic" values are R = (l-3) x 10"'^ 
cm^ s~' and / = 10"** photons cm~^ s~' Hz~'. In Table 7 we de- 
fine four models grids for comparison with the observed abun- 
dance and excitation patterns. These grids are chosen to have 
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the full range of Tidn, n, and d but vary in their assumed incident 
radiation and H2 formation rate constants. 

In Figure | we display the LMC, SMC, and S77 samples to- 
gether with the model clouds, in terms of molecular fraction 
correlated with total H content. First, we compare the data 
to model grid A, with typical Galactic conditions (upper left 
panel). These models match the Galactic points reasonably 
well, but they do not coincide with the LMC and SMC sam- 
ples at all. The Milky Way agreement is a good starting point 
for understanding the different conditions in the low-metallicity 
Clouds. 

In Figure ^ (upper right panel), we plot model grid B, with a 
typical Galactic radiation field and a low formation rate con- 
stant, R = 3 X 10"'^ cm-' s"'. This value is 3-10 times be- 
low the Galactic values inferred from Copernicus data, R = 
(1 -3) X 10"'^ cm^ s"' (Jura 1974). These model clouds match 
some LMC points, but in general they are overabundant in H2, 
particularly at the column densities N(H 1) of the SMC sight 
lines. In the lower left panel, we show the effect of enhancing 
the radiation field incident on the model clouds to / = 10"^- 10"*' 
cm"^ s"' Hz"', 10-100 times the Galactic mean (Jura 1974), 
while fixing R at the Galactic value (model grid C). This ad- 
justment again produces the correct change in the pattern, re- 
ducing the molecular fraction and coinciding with some LMC 
points. Finally, the lower right panel displays model grid D, 
raising the radiation field by factors of 10 - 100 and lowering 
R to 1/3 - 1/10 the Galactic value. Only in this extreme case 
are we able to produce an abundance pattern that resembles the 
SMC and LMC samples. The observed clouds have relatively 
high N(H I), but lower average /h2 than the S77 sample. A 
combination of reduced R and enhanced /, though not unique, 
can reproduce the observed patterns. 

We have made no attempt to derive individual model param- 
eters for each of the sight lines in the sample. Such solutions 
are not unique in the absence of independent constraints on the 
cloud density and kinetic temperature (Browning, Tumlinson, 
& Shull 2001, in preparation). Indeed, such sight-line model- 
ing is unnecessary for a large sample that can be compared to 
similar ensembles of model clouds. The non-uniqueness prob- 
lem and the large variation in the observed conditions implied 
by the LMC and SMC abundance distributions hinder us from 
deriving unique quantitative constraints on j3 and R for the sight 
lines. However, we have shown that realistic grids of model 
clouds that accurately reproduce well-known Galactic patterns 
can match the LMC/SMC molecular abundance patterns seen 
in the FUSE samples. These models assume high incident radi- 
ation fields and a low grain formation rate coefficient. 

3.3. H2 Excitation 

We cannot determine from the abundance information alone 
whether the data favor suppressed molecule formation in low- 
metallicity gas or enhanced photo-dissociation by UV radiation 
as an explanation for the reduced (/) . The most likely scenario 
is a combination of these two effects. To separately diagnose 
the effects of the competing formation/destruction processes, 
we examine the rotational excitation of the H2. 

Following the absorption of a FUV photon and excitation to 
an upper electronic state, the H2 molecule fluoresces to the 
ground electronic state and cascades down through the rota- 
tional and vibrational lines. The net effect of these repeated 
excitations and cascades is a redistribution of the molecules into 
the excited rotational states of the ground electronic and ground 
vibrational state. We observe this distribution of molecules di- 



rectly in our data, and derive column densities N(J) in each ro- 
tational level of the ground vibrational state. 
The rotational temperature Tqi is given by: 



mgilgQ)N{Q)lN{l)V 



(5) 



where the column densities are those given in Tables 5 and 6, 
the ratio gi/gQof statistical weights of J = 1 and J = is 9, and 
AZioi /k=\l\K. The quantity Tq\ is thought to trace the kinetic 
temperature of the molecular gas at high densities where colli- 
sions dominate the level populations (such as the translucent 
cloud regime, see Snow et al. 2000). In lower-column diffuse 
clouds, the relationship between Tqi and Tkin is not known ex- 
actly. As a simple test for widely divergent physical conditions 
in the Magellanic Clouds relative to the Galaxy, we derive the 
rotational temperature T^i for our targets (Tables 3 and 4). In 
Figure [l^ (top left panel) we plot total H2 column density ver- 
sus T()\. A clear pattern emerges, with widely varying Tqi at 
lower N(H2), where clouds are not expected to self-shield and 
densities may be too low for collisions with H""" or H3 to estab- 
lish a thermal ortho-para ratio. For all sight lines with N(H2) 
> 10'^-^ cm-2, we find {Tm) = 82±21 K (the standard devi- 
ation is calculated excluding the outlier AV 47 in the SMC). 
For all sight lines, we find {Tq\) = 1 15 K. Thus, for likely self- 
shielded clouds, the rotational temperatures lie near the Galac- 
tic average, 77 ± 17 K (S77), providing evidence that the dis- 
tribution of kinetic temperatures is similar to that in Galactic 
diffuse molecular gas. While Tq\ may closely track the kinetic 
temperature, the measured Doppler b parameters indicate sig- 
nificant non-thermal motions and do not correlate with Tqi. 

To further explore the rotational excitation of H2 in the sam- 
ples, we use the column density ratios N(3)/N(l), N(4)/N(2), 
and N(5)/N(3), in order of decreasing sensitivity to coUisional 
excitation and increasing sensitivity to the radiative cascade. 
We plot these ratios for the LMC, SMC, and Galactic compari- 
son samples in Figure 10 The trend toward lower excitation at 
higher N(H2) is apparent in all the ratios. In N(3)/N(l), there 
is no clear distinction between the three samples, perhaps in- 
dicating the importance of collisions in determining this ratio. 
However, it appears that both the LMC and SMC samples are 
enhanced in N(4)/N(2) and N(5)/N(3), perhaps indicating the 
presence of a FUV radiation field above the Galactic mean. The 
large error bars on these ratios are a direct result of the uncer- 
tainties in the Doppler b parameter. 

In Figure |ll] we again plot the column density ratios 
N(4)/N(2) and N(5)/N(3) compared with the model grids. In 
this comparison, the difference between the Galactic and Mag- 
ellanic samples is striking. With the exception of six outliers, 
all the Galactic stars lie within the region described by model 
grids A and B. In contrast, the majority of the LMC and SMC 
points have N(4)/N(2) and N(5)/N(3) higher than the model grid 
with Galactic conditions. 

We consider four hypotheses to explain the unusually high 
excitation ratios seen in the Clouds: (1) low kinetic tempera- 
tures elevate the ratios by differentially enhancing J = over 
J = 1, which are then redistributed to J = 4 and 5, respec- 
tively; (2) the upper levels are preferentially populated in newly 
formed molecules ("formation pumping"); (3) high incident 
FUV radiation pumps the upper levels via a radiative cascade; 
(4) a low formation rate of H2 on grain surfaces allows FUV 
radiation to penetrate to larger N(H2), elevating the ratios in 
the cloud cores. The results of our modeling follow these four 
mechanisms. 
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(1) Low kinetic temperature elevates the excitation ratios in 
model clouds at fixed density, size, and incident radiation. The 
J = 4 and 5 levels are more efficiently populated than J = 2 and 
3 by absorptions from J = and 1, followed by radiative cas- 
cade. Thus, any process that boosts N(0) and N(l) increases 
N(4)/N(2) and N(5)/N(3). The right envelope of the models 
pictured in the upper left panel of Figure |ll] is defined by r^in = 
10 K. Low r^in by itself is insufficient to explain the enhanced 
ratios. 

(2) The idea that H2 molecules are formed on the surfaces of 
interstellar dust grains has widespread theoretical and observa- 
tional support. However, the rotational and vibrational excita- 
tion of the H2 molecule as it exits the grain is unknown. Be- 
cause of the high Einstein A values of the levels involved in the 
initial fluorescence, all memory of the formation distribution is 
quickly erased in irradiated ensembles of molecules. The cur- 
rent numerical code assumes that all molecules are formed in a 
3:1 ortho-para ratio in the highest vibrational level (v = 14) of 
the ground electronic state. This assumption was also made by 
the models of van Dishoeck & Black (1986). However, because 
in equilibrium there is only one molecule formed for every 10 
excitations by a UV photon (fdiss ^ 0.1), the formation distri- 
bution is quickly erased. Indeed, test models produced by the 
current code over the range of radiation inputs / and assuming 
that half of all H2 is formed into J = 4 and half into J = 5 do 
not differ appreciably from models made with this general as- 
sumption. Thus, formation pumping is an unlikely cause of the 
elevated excitation. 

(3) As end points to the radiative cascade from the upper 
electronic state, J = 4 and J = 5 levels are favored over J = 2 and 
J = 3. Thus, more frequent photoabsorptions result in higher 
N(4)/N(2) and N(5)/N(3) in model clouds. In Figure [u] (lower 
panel) we show the complete grid of models, which range up 
to /= 1 .0 X 10"'' photons cm"^ s"' Hz"' , 100 times the Galactic 
mean (Jura 1974). The points with Tkin = 10 K and this high ra- 
diation field are still not consistent with the observed excitation 
ratios. 

(4) A low grain formation rate of H2 can elevate the exci- 
tation ratios in model clouds. A smaller formation rate allows 
radiation to penetrate to higher N(H2) in deeper parts of the 
cloud, similar to the effects of increasing the incident radiation 
field. However, models with R = 3 x 10"'** cm^ s"', one-tenth 
the typical Galactic value, do not reproduce the observed ratios. 

Finally, we attempt to match the model clouds with the same 
set of parameters that successfully match the H2 abundance pat- 
terns. The lower panels of Figure 111] show the grid of 1890 
models with R = 3x 10"'^ and / = 10"^ - 10"*". The combina- 
tion of low R and high / produces the best match to the ob- 
served N(4)/N(2) ratios and the N(5)/N(3) ratios at low N(H2). 
Based on this concordance and the match to the molecular frac- 
tion results above, we conclude that, on average, diffuse H2 in 
the Magellanic Clouds exists in conditions different from the 
Galaxy. These sight line models are not unique, but both a high 
radiation field and a low formation rate coefficient are neces- 
sary to explain the observed abundance and excitation patterns 
in grids of model clouds that accurately reproduce the Galactic 
distributions. 

The model grids with high / and low R provide the best match 
to the observed abundance and excitation of H2, but there are 
still some LMC and SMC points with discrepant excitation ra- 
tios. If more than one distinct H2 cloud exists in the sight line, 
their column densities can sum and give anomalous column 
density ratios that are not reproducible by one-sided slab mod- 



els. Large N(4)/N(2) can occur, for instance, if two clouds are 
present: one cloud that contributes the majority of the H2 col- 
umn density but no N(4), and one that contributes little N(H2), 
little N(2), but N(5) ^ lO'"*-^ cm"^. This combination of clouds 
can appear in Figure 11 to have log[N(4)/N(2)] = -1.0, while 
no single cloud in the model grid can achieve this high ratio. 
However, it is unlikely that concatenation of multiple clouds 
will occur more frequently in the LMC and SMC than in the 
Galaxy and thus explain their anomalous column density ratios. 

As we noted in the discussion of the H2 abundance, indi- 
vidual sight line models are non-unique. However, taken to- 
gether, the reduced abundance and enhanced rotational excita- 
tion of the LMC and SMC H2 is strong evidence for a shift in 
the formation/destruction balance in the low-metallicity gas of 
the Clouds. Grain formation rates 1/3 to 1/10 times the Galactic 
value and UV photo-absorption rates 10 to 100 times the Galac- 
tic mean value are necessary to explain the observed patterns of 
H2 abundance and rotational excitation. 

3.4. Dust and Gas Correlations 

From the correlations of H2 and H I in the samples, it is clear 
that large increases in the fractional gas content of the Mag- 
ellanic Clouds relative to the Galaxy have not led to a corre- 
sponding increase in the amount of molecular gas in the diffuse 
ISM. This fact, and the lower dust-to-gas ratios in the Clouds, 
confirm that interstellar dust plays a large role in the production 
of diffuse H2. Thus, correlations with the dust content of the 
clouds we probe may reveal such a link between dust and H2. 
We use dust emission to test for a correlation between H2 and 



dust abundance . Figures 1 2 and 1 3 show IRAS 1 00 //m maps of 



the LMC and SMC, overlaid with surface brightness contours 
at 10, 100, and 1000 MJy sr"'. Sight fines in which we detect 
H2 are marked with +, while sight lines for which we do not de- 
tect H2 are marked with x . For the LMC detections, the IRAS 
100 /im surface brightness ranges from 6 - 2240 MJy sr"' with 
a mean of 277 MJy sr"' . The LMC non-detections have 7 - 146 
MJy sr"' with a mean of 51 MJy sr"'. For the SMC, the H2 
detections have IRAS brightnesses of 3 - 110 MJy sr"' with a 
mean of 39 MJy sr"'. The two SMC non-detections have FIR 
surface brightness ctfir = 23 and 40 MJy sr"' . Thus, the ob- 
served far-infrared flux at a given position is not a predictor of 
detectable H2, except in regions of strong dust emission such 
as 30 Doradus. However, we note that the 1.5' IRAS pixel size 
subtends ^ 22 pc at the distance of the LMC and 26 pc at the 
SMC. Shull et al. (2000) inferred linear sizes t/ ~ 2-3 pc for 
three Galactic H2-bearing clouds. The different molecular ex- 
citation in the HD 5980 and AV 232 sight lines, separated by 
17 pc at the SMC, argues for clouds sizes smaller than the size 
scale probed by the IRAS pixel. These two techniques may 
sample different distributions of gas in the ISM, and perhaps 
should not be expected to correlate. 

3.5. Diffuse H2 Mass of the LMC and SMC 

Over the past decade, the Magellanic Clouds have been stud- 
ied in increasing detail in mm-wave CO emission. In all cases, 
the results point to CO emission weaker by factors of 3-5 than 
expected if Milky Way GMCs were observed at a distance of 
50-60 kpc. Israel et al. (1982, 1986) studied the LMC and 
SMC in '^CO (2-1) emission, using the ESO 3.6m telescope 
at La Silla, Chile with a 2' FWHM beam at 230 GHz (1.3 mm). 
They found relatively weak CO emission, with peak antenna 
temperatures about 25% of the expected signal from typical 
Galactic molecular clouds at Magellanic Cloud distances. They 
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Table 7 
H2 Model Cloud Grids 



Label 




7(10-**) 


Description 




^ -1 
cm s 


ph cm"^ s"' Hz"' 




A 


1-3 


1 -4 


Galactic conditions 


B 


0.3 


1 -4 




C 


1-3 


10-100 




D 


0.3 


10-100 


LMC, SMC 



also found that the CO emission is less widespread in the SMC, 
compared to either the Galaxy or the LMC. More recent sur- 
veys with the ESO/SEST telescope at 45" resolution (Israel et 
al. 1993; Rubio, Lequeux, & Boulanger 1993) and the NAN- 
TEN 4 m telescope at 2.6' resolution (Fukui et al. 1999) find 
similarly weak '^CO (1-0) emission in the LMC and SMC. 

The preliminary estimate for the calibration between '^CO 
line intensity and total hydrogen column density in the SMC 
(Israel et al. 1993) is approximately four times larger than 
the canonical Galactic value, Xg^i = [N(H I)h-2N(H2)]//co ~ 
2 X 10^° cm-2 (K km s"')-'. For the LMC, Fukui et al. (1999) 
find a ratio Xlmc = 9 x 10^" cm"^ (K km s"')"', about half the 
value derived by Israel (1997) using IRAS emission. The NAN- 
TEN survey estimates a total H2 cloud mass of (4-7) x 10^ Mq 
for clouds with N(H2) > 1 x 10^' cm"^, corresponding to 5 - 
10% of the atomic mass. 

Most of these studies (Israel et al. 1986; Rubio et al. 1993) 
conclude that the weak CO emission in the Magellanic Clouds 
does not arise from a deficit of molecular gas, but rather from 
the low metallicity and high UV radiation fields. These fac- 
tors result in low C and O abundances and high CO photodis- 
sociation rates, which reduce the CO/H2 ratio from the stan- 
dard Galactic values. These suggestions have been sustained 
by plane-parallel models of molecule formation and radiative 
transfer (Maloney & Black 1988; Lequeux et al. 1994). In 
plane-parallel models of clouds with N(H2) = 10^^ cm"^, Mal- 
oney & Black (1988) find that the GMCs with low (SMC) 
metallicity and correspondingly reduced dust/gas ratios should 
have much smaller CO column densities and CO-emitting sur- 
face areas. They also find that the H2 number densities are not 
as sensitive to the lower metallicities. Consequently, the LMC 
and SMC molecular clouds are expected to show low CO/H2 
ratios. In this section, we use the FUSE data on H2 to directly 
address the issue of the H2 mass fractions in the LMC and SMC. 

Using the observed distribution of H2 and the high-resolution 
H I emission maps of Stanimirovic et al. (1999) and Staveley- 
Smith et al. (2001), we can estimate the total mass of H2 in the 
diffuse regions of the LMC and SMC. To assign an H2 column 
density to each observed position in the H I column density 
map, we use a plot of N(H2) versus N(H I) for the FUSE targets 
to produce 5 bins of ANhi = 1 x 10^' cm"^. We assign a molec- 
ular fraction fui to each point in the bin, with fn2 randomly 
selected from the set of FUSE targets with N(H I) in that bin. 
This technique is used to produce a synthetic H2 map. A direct 
integral over this map produces a total diffuse H2 mass for the 
Clouds. 

For the LMC, we obtain a total atomic hydrogen mass 



4.8 X 10 Mq. From the H I map and the observed FUSE dis- 
tribution of H2, we obtain a molecular mass N(H2) = 8 x 10^ 
Mq. Thus, we find that ^ 2% of the hydrogen mass in the 
LMC resides in diffuse H2. This is consistent with the cumu- 
lative molecular fraction for the LMC, (fm) = 0.012. For the 
SMC, this technique determines a total atomic hydrogen mass 
4 X 10** Mq. We find a molecular mass N(H2) = 2 x 10* M©. 
Thus, 0.5% of the hydrogen mass of the SMC lies in the diffuse 
H2. 

Our inferred masses represent less than 1% of the diffuse H I 
measured in the Clouds. Israel (1997) used IRAS 100 fim maps 
and detected CO clouds (Cohen et al. 1988) in the LMC to pre- 
dict the surface density of H2 in the LMC as a function of CO 
surface brightness in radio measurements. That study found 
total molecular masses of (1.0 ±0.3) x 10** M© for the LMC 
and (0.75 ±0.25) x 10^ Mq for the SMC. These masses cor- 
respond to ~ 20% mass fraction in the Clouds, 10 times that 
found by the FUSE survey. As noted above, the NANTEN sur- 
vey (Fukui et al. 1999) found values about half this amount for 
the LMC. We compare the detected column densities of H2 to 
these predictions in Figure |l4 . Although we do not probe di- 
rectly the CO clouds on which the Israel (1997) and Fukui et 
al. (1999) results are indirectly based, the FUSE target stars do 
coincide with regions of similar IRAS 100 /im surface bright- 
ness. We find that total H content in our FUSE sight lines is 
well correlated with the FIR surface brightness, but in the range 
of overlap, our detected H2 column densities lie 1-5 orders of 
magnitude below that predicted by the IRAS flux. Only one 
FUSE sight line coincides with the predicted surface density, 
indicating that the total molecular mass in the LMC may be far 
lower than that predicted by the total dust and/or CO content. 
Based on Figure |l^ we conclude that the relationships between 
N(CO), N(H I), and FIR surface brightness that predict high 
N(H2) do not apply in diffuse regions away from the detected 
CO clouds. These regions, which we probe directly with the 
FUSE sight Unes, possess a wide range of molecular fraction 
but are generally below the value fm ^0.1 required to produce 
M(H2) / M(H I) = 0.2. 

We derived M(H2) for the Clouds by assigning observed 
molecular fractions to the 21 cm emission maps according to 
the observed H I column density. In both the SMC and LMC 
maps, a small fraction of the beam positions have measured 
N(H I) larger than any H I column probed in this survey. In 
the LMC, the maximum N(H I) exceeds the highest value in 
our survey by 2x, and positions with N(H I) > 5.0 x 10^' cm"^ 
comprise 40% of the H I by mass. In the SMC, N(H I) exceeds 
the survey maximum by 30%, but these positions comprise only 
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0.4% of the H I by mass. Thus, in the SMC, this survey has 
sampled the distribution of observed N(H I) rather completely, 
while in the LMC the upper end of the distribution has not been 
sampled. 

The discrepancy between H2 masses determined from FUSE 
and those predicted from the 100 ^m flux can be resolved if 
the H I positions with N(H I) higher than our sample are highly 
molecular. However, the molecular fractions in these points 
must be extremely high to achieve M(H2) ^10^ Mq. For the 
SMC, those points with N(H I) > 8.0 x 10^' cm"^ must have 
/h2 = 0.87 to give M(H2) = 7.5 x 10^ Mq. For the LMC, where 
N(H I) ranges to 2x that in our sample, /h2 = 0.70 is required. 
These molecular fractions are substantially higher than the fm 
^0.1 seen in high column-density Galactic molecular clouds 
and are somewhat like the /h2 ^ 0.7 seen in Galactic translu- 
cent clouds (Snow et al. 2000; Rachford et al. 2001) or GMCs. 
Thus the molecular mass discrepancies are resolved since our 
survey does not sample the upper end of the molecular cloud 
regime. A large fraction of H2 could reside in cold, dense gas 
that appears in neither the H I emission map nor in the FUSE 
data, which probes to E'(B-V) = 0.20 only. This question will 
be addressed by FUSE Cycle 2 observations of seven sight lines 
to the SMC and LMC wifli E(B-V) > 0.20. 

4. DISCUSSION AND CONCLUSIONS 

We have surveyed the abundance, excitation, and properties 
of diffuse H2 in the diffuse interstellar medium of the Small 
and Large Magellanic Clouds. The central results of this sur- 
vey are the column density measurements themselves (Tables 
5 and 6) and the statistics of H2 detections and non-detections. 
In addition to these basic findings, we have also correlated the 
detected H2 with the H I and dust along these sight lines, and 
derived rotational excitation parameters of the H2. 

Our major conclusions are: 

• In a sample of 70 sight lines to hot stars in the LMC and 
SMC, we find H2 towards 23 of 44 LMC stars and to- 
wards 24 of 26 SMC stars. 

• The overall abundance of diffuse H2 in the ISM of the 
Magellanic Clouds is lower than in the Milky Way. We 
find molecular fractions (/h2 



0.010![!;JJ[;^ for the SMC 



and {fm) = 0.012!;j;™f for the LMC, compared with 
ifm) = 0.095 for the Galactic disk over a similar range of 
reddening. The main source of error is systematic uncer- 
tainty in the measurement of N(H I). This result can be 
interpreted as the effect of suppressed H2 formation on 



dust grains and enhanced H2 destruction by FUV pho- 
tons. 

• The reduced molecular fraction and enhanced rotational 
excitation in the Clouds are consistent with models hav- 
ing a low H2 formation rate coefficient and a high UV 
photo-absorption rate, relative to Galactic conditions. 
Using grids of models clouds, we find that ^ 3 x 10"'^ 
cm^ s-\ 1/3 - 1/10 the Galactic rate, and / - 10"^- lO"*" 
photons cm"^ s"' Hz"', 10 to 100 times the Galactic 
mean, are necessary to replicate the LMC and SMC re- 
sults. 

• The mean molecular fraction of the detected diffuse gas 
implies a total diffuse H2 mass of M(H2) = 8x10* Mq 
for the LMC and M(H2) = 2 x 10* Mq for the SMC. 
These masses are <^ 2% of the H I masses of the Clouds, 
implying small overall molecular content, high star for- 
mation efficiency, and/or substantial molecular mass in 
cold, dense clouds unseen by our survey. 

• We find no correlation between observed H2 and dust 
properties below a critical value of E(B-V) and/or a crit- 
ical 100 /im IRAS flux. However, high extinction and/or 
strong dust emission are predictors of high H2 abun- 
dance. 

This survey is one of four major surveys of interstellar H2 
that are being conducted by the Colorado group with FUSE. 
The Galactic H2 survey consists of over 100 sight lines to hot 
stars in the disk, with \b\ < 10°. The Galactic halo survey con- 
sists of '^lOO sight lines to the Magellanic Clouds and extra- 
galactic AGN and QSOs. Finally, the FUSE translucent cloud 
survey consists of 35 stars in the Galactic disk, with Ay ^ 2.0. 
Taken together, these surveys and the detailed followups should 
form the basis of a more complete understanding of molecule 
formation in the diffuse ISM than has been possible to date. 

Lister Staveley-Smith graciously provided the LMC H I 
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Space Flight Center We have also made extensive use of the 
SIMBAD database, operated at CDS, Strasbourg, France. This 
work was supported in part by the Colorado astrophysical the- 
ory program through NASA grants NAG5-4063 and NAG5- 
7262. 
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Table 1 

FUSE H2 Large Magellanic Cloud Target List 



IN dine 




UGL' yjz^yjyjyj) 


opeciidi lypc 


Jvci. 


^..'Dservaiiun 




f 






h m s 


d m s 






LL) 


(used/total) 


(ks) 


(per pixel)'' 


Sk -67 05 


4 50 18.96 


-67 39 37.90 


09.7 lb 


1 


P1030703 


1/1 


8.0 


10 


Sk -67 14 


4 54 31.92 


-67 15 24.90 


B1.5Ia 


4 


PI 174201,3 


3/5 


10.5 


7 


HD 32109 


455 31.50 


-67 30 1.00 


WN4b 


5 


PI 174401-2 


5/5 


16.7 


8 


Sk -65 21 


5 1 22.33 


-65 41 48.10 


09.7 lab 


8 


P1030901-4 


4/4 


18.2 


14 


Sk -65 22 


5 1 23.14 


-65 52 33.50 


06 Iaf+ 


1 


P1031002 


11/11 


27.2 


9 


HD 33133 


5 3 9.00 


-66 40 57.76 


WN8h 


5 


PI 174501 


4/4 


4.6 


7 


Sk -69 59 


5 3 12.79 


-69 1 37.20 


BO 


2 


P1031103 


4/4 


25.7 


9 


Sk -70 60 


5 4 40.94 


-70 15 34.50 


04-5 V:n 


9 


PI 172001 


2/2 


7.9 


6 


Sk -70 69 


5 5 18.73 


-70 25 49.80 


05 V 


8 


P1172101 


2/2 


6.1 


6 


Sk -68 41 


5 5 27.20 


-68 10 2.70 


B0.5 la 


4 


PI 174102 


1/1 


6.5 


6 


Sk -68 52 


5 7 20.60 


-68 32 9.60 


BO la 


1 


PI 174001 


2/5 


8.1 


7 


Sk -67 69 


5 14 20.16 


-67 8 3.50 


04 Ill(f) 


7 


P1171703 


4/4 


7.8 


7 


Sk-69 104 


5 18 59.57 


-69 12 54.70 


06 Ib(f) 


1 


PI 172401 


1/1 


4.1 


9 


Sk -67 76 


5 20 5.81 


-67 21 8.90 


BO 


2 


P1031201 


3/5 


19.2 


9 


Sk -69 124 


5 25 18.37 


-69 3 11.10 


09 lb 


10 


PI 173601-2 


3/3 


12.7 


9 


Sk -67 101 


5 25 56.36 


-67 30 28.70 


08 II((f)) 


9 


PI 173402 


2/2 


6.2 


6 


Sk -67 104 


5 26 4.15 


-67 29 56.50 


WC4(H-0?)+08I: 


3 


P1031302 


3/3 


5.1 


9 


Sk -68 80 


5 26 30.43 


-68 50 26.60 


WC4H-06V-III 


3 


P1031402 


4/4 


9.7 


9 


BI 170 


5 26 47.79 


-69 6 11.70 


09. 5 lb 


9 


PI 173701 


1/1 


4.3 


7 


Sk -67 111 


5 26 48.00 


-67 29 33.00 


06: lafpe 


9 


PI 173001 


2/2 


8.0 


8 


BI 173 


5 27 10.08 


-69 7 56.20 


08 II: 


9 


PI 173201-2 


4/4 


11.4 


8 


Sk -70 91 


5 27 33.74 


-70 36 48.30 


06.5 V 


10 


PI 172501 


1/1 


5.5 


7 


Sk -66 100 


5 27 45.59 


-66 55 15.00 


06 n(f) 


8 


PI 172303 


2/2 


7.1 


7 


HD 269582 


5 27 52.75 


-68 59 8.60 


WNlOh 


11 


PI 174701 


1/1 


4.6 


7 


HD 37026 


5 30 12.22 


-67 26 8.40 


WC4 


12 


P1175101 


3/3 


8.1 


6 


Sk -71 45 


5 31 15.55 


-71 4 8.90 


04-5 Ill(f) 


1 


P1031501-4 


4/4 


18.9 


15 


HD 269687 


5 31 25.61 


-69 5 38.40 


WNllh 


11 


PI 174801 


1/1 


3.8 


6 


Sk -67 166 


5 31 42.00 


-67 38 6.00 


04If+ 


1 


A1330100 


84 


220. 


65 


Sk -67 169 


5 31 51.68 


-67 2 22.30 


Bl la 


4 


P1031603 


4/7 


40.5 


8 


Sk -67 167 


5 31 51.98 


-67 39 41.10 


04 Inf+ 


7 


PI 171902 


1/1 


2.8 


10 


Sk -67 191 


5 33 34.12 


-67 30 19.60 


08 V 


10 


P1173102 


1/1 


6.6 


6 


BI 208 


5 33 57.45 


-67 24 20.00 


07 Vz 


9 


PI 172702-5 


4/4 


15.6 


9 


HD 37680 


5 34 19.39 


-69 45 10.00 


WC4 


12 


PI 175001 


3/3 


7.0 


5 


HD 269810 


5 35 13.92 


-67 33 27.00 


02 III(f*) 


6 


PI 171603 


5/5 


8.2 


9 


BI 229 


5 35 32.20 


-66 2 37.60 


07 V-III 


9 


PI 172801 


2/2 


5.6 


7 


Sk -66 169 


5 36 54.50 


-66 38 25.00 


09.7 la-H 


4 


P1173801 


3/3 


4.8 


8 


Sk -66 172 


5 37 5.56 


-66 21 35.70 


02in(f*)H-0B 


6 


PI 172201 


1/1 


3.6 


6 


Sk -68 135 


5 37 48.60 


-68 55 8.00 


ON9.7 Ia+ 


1 


PI 173901 


3/3 


6.8 


5 


Mk42 


5 38 42.10 


-69 5 54.70 


03 If*AVN6-A 


14 


PI 171802-4 


5/5 


15.4 


15 


Sk -69 243 


5 38 42.57 


-69 6 3.20 


WN+OB composite 


13 


P103 1705-6 


2/2 


9.5 


12 


Sk -69 246 


5 38 53.50 


-69 2 0.70 


WN6h 


5 


P1031802 


4/4 


22.1 


9 


HD 269927 


5 38 58.25 


-69 29 19.10 


WN9h 


5 


PI 174601 


2/2 


7.2 


8 


BI 272 


5 44 23.18 


-67 14 29.30 


07: III-II: 


9 


PI 172902 


1/1 


6.6 


7 


Sk -70 115 


5 48 49.76 


-70 3 57.50 


06.5 laf 


9 


PI 172601 


2/2 


4.8 


7 



^Spectral type references: (1) Walborn (1977); (2) Rousseau et al. (1978); (3) Moffat et al. (1990); (4) Fitzpatrick (1988); 
(5) Smith et al. (1996); (6) Walborn et al. (2001b); (7) Garmany & Walborn (1987); (8) Walborn et al. (1995); (9) Walborn et 
al. (2001a) (10) Conti et al. (1986); (11) Crowther & Smith (1997); (12) Torres-Dodgen & Massey (1988); (13) Breysacher et 

al. (1999); (14) Walborn & Blades 1997. 

''Signal-to-noise ratio per pixel in the 1040-1050 A range. Signal-to-noise ratio per resolution element varies with the spectral 
resolution, which is not fixed in our survey. 
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Table 2 

FUSE H2 Small Magellanic Cloud Target List 



Name 


RA (J2000) 


Dec (J2000) 


Spectral Type 


Ref.'' 


Observation 


Exposures 


^exp 


S/N 




h m s 


d m s 






ID 


(used/total) 


(ks) 


(per pixel)'' 


AV 14 


46 32.66 


-73 6 5.60 


03-4 V 


1 


PI 175301 


2/2 


6.8 


5 


AV 15 


46 42.19 


-73 24 54.70 


06.5 11(f) 


2 


P1150101 


12/12 


14.6 


7 


AV26 


47 50.07 


-73 8 20.70 


07 III 


1 


PI 176001 


1/1 


4.0 


4 


AV47 


48 51.35 


-73 25 57.60 


08 m((f)) 


2 


PI 150202 


5/5 


16.3 


8 


AV69 


50 17.40 


-72 53 29.90 


OC7.5 III((f)) 


2 


PI 150303 


13/13 


17.6 


9 


AV 75 


50 32.50 


-72 52 36.20 


05 III(f+) 


2 


PI 150404 


4/4 


14.4 


7 


AV83 


50 52.01 


-72 42 14.50 


07Iaf+ 


2 


PI 176201 


1/1 


4.0 


5 


AV95 


51 21.54 


-72 44 12.90 


07 mm 


2 


PI 150505 


11/11 


14.3 


7 


AV207 


58 33.19 


-71 55 46.50 


07 V 


3 


PI 175901 


3/3 


4.1 


3 


NGC 346-6 


58 57.74 


-72 10 33.60 


04V((f)) 


6 


PI 175601 


4/4 


6.4 


8 


NGC 346-4 


59 0.39 


-72 10 37.90 


05-6 V 


5 


PI 175701 


2/2 


5.2 


8 


NGC 346-3 


59 1.09 


-72 10 28.20 


02 in(f*) 


12 


PI 175201 


2/2 


5.1 


8 


HD 5980 


59 26.57 


-72 9 53.90 


WNvar 


4 


X0240202 


8/8 


3.2 


9 


AV232 


59 32.19 


-72 10 46.20 


07 Iaf+ 


8 


X0200201 


5/5 


10.3 


8 


Sk 82 


59 45.72 


-72 44 56.10 


BO law 


9 


P1030301 


10/10 


16.2 


8 


AV238 


59 55.61 


-72 13 37.70 


09.5 III 


2 


PI 176601 


2/2 


11.1 


6 


AV242 


1 6.84 


-72 13 57.00 


Bl la 


7 


PI 176901 


2/2 


5.1 


6 


AV264 


1 1 7.72 


-71 59 58.60 


Bl la 


7 


PI 177001 


3/3 


4.5 


5 


AV321 


1 2 57.04 


-72 8 9.30 


BO IIIww 


1 


PI 150606 


6/6 


16.9 


7 


AV327 


1 3 10.58 


-72 2 13.80 


09.5 Il-Ibw 


2 


PI 176401 


4/4 


4.7 


6 


Sk 108 


1 3 25.24 


-72 6 43.30 


WN3 + 06.5(n) 


8 


X0150102 


8/8 


13.3 


4 


AV372 


1 4 55.73 


-72 46 47.70 


09.5 labw 


10 


PI 176501 


1/1 


4.4 


5 


AV378 


1 5 9.44 


-72 5 35.00 


08 V 


1 


PI 150707 


8/8 


14.7 


7 


AV388 


1 5 39.62 


-72 29 26.80 


04 V 


1 


PI 175401 


2/2 


5.6 


6 


Sk 159 


1 15 58.84 


-73 21 24.10 


B0.5 law 


9 


P1030501 


2/2 


6.9 


7 


Sk 188 


1 31 4.23 


-73 25 2.20 


W03-H04V 


11 


P1030601-4 


8/8 


23.1 


15 



^Spectral type references: (1) Garmany et al. (1987); (2) Walborn et al. (2000); (3) Crampton & Greasley (1982); (4) 
Koenigsbergeret al. (1994); (5) Walborn & Blades (1986); (6) Walborn et al. (1995); (7) Lennon (1997); (8) Walborn (1977); 
(9) Walborn (1983); (10) Walborn et al. (2001a); (11) Crowther et al. (1998); (12) Walborn et al. (2001b). 

''Signal-to-noise ratio per pixel in 1040-1050 A range. Signal-to-noise ratio per resolution element varies with spectral 
resolution, which is not fixed in our survey. 
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Table 3 

FUSE H2 Large Magellanic Cloud Summary 



Name E(B-V)^ vlsr'' N(H2) Jqi ^oopp N(HI)^ 
(kms-i) (cm-2) (K) (kms"!) (lO^o cm-^) 



Sk -67 05 


0.16 


256 ±7 


Sk -67 14 


0.08 


268 ±3 


HD 32109 


0.05 


284 ± 13 


Sk -65 21 


0.12 


222 ± 11 


Sk -65 22 


0.15 


269 ±4 


HD 33133 


0.08 


299 ± 14 


Sk -69 59 


0.12 




Sk -70 60 


0.135 




Sk -70 69 


0.09 




Sk -68 41 


0.105 




Sk -68 52 


0.17 


232 ± 12 


Sk -67 69 


0.15 


278 ± 11 


Sk -67 104 


0.13 




Sk -67 76 


0.11 


247 ±8 


Sk -69 124 


0.125 




Sk-67 101 


0.115 




Sk -69 104 


0.105 




Sk -68 80 


0.10 




BI 170 


0.125 




Sk-67 111 


0.115 




BI 173 


0.16 


217± 11 


Sk -70 91 


0.10 


281 ± 8 


Sk-66 100 


0.11 




HD 269582 


0.09 




HD 37026 


0.11 


246 ±9 


Sk -71 45 


0.18 


228 ±9 


HD 269687 


0.10 




Sk -67 166 


0.10 


209 ±2 


Sk -67 169 


0.07 




Sk-67 167 


0.12 




Sk-67 191 


0.11 




BI 208 


0.03 


265 ± 11 


HD 37680 


0.15 


222 ± 11 


HD 269810 


0.10 




BI 229 


0.135 




Sk -66 169 


0.14 




Sk -66 172 


0.19 


267 ±9 


Sk-68 135 


0.27 


250 ±5 


Sk -69 246 


0.10 


261 ±9 


HD 269927 


0.19 


257 ±7 


BI 272 


0.14 




Sk-70 115 


0.21 


202 ± 4 



19.46^0:05 


57 


6.2!}:1 


7.7 


15.01!°:?^ 


270 


linear 


8.7 


18.67!°;! 


72 
90 


2.3!?:? 


5.6 
5.4 


14 93+0.19 


117 


linear 


10.3 


1 7 cn+0.26 


30 


3 5+2.6 


8.6 


<14.10 






10.2 


< 14.46 






3.5 


<14.56 






3.3 


< 14.30 






10.8 




60 


4.8!?| 


12.8 


15.47!j:5i 


115 


6.9!g 


12.5 


<14.15 






13.6 


14.97!jj-^^ 


120 


linear 


7.9 


< 14.26 






8.3 


< 14.26 






13.4 


< 14.40 






13.8 


<14.10 






7.1 


< 14.26 






7.2 


< 14.86 






13.9 


15.64!«:f3 


117 


c Q+2.2 
-^.^-1.7 


6.2 


i4.8i!«:| 


156 


linear 


6.2 


< 14.40 






3.2 


< 14.30 






6.8 


1 ^ A^+l 94 
1J.'+J_0;25 


91 


5.4!^:^ 


7.3 


1 09 

i8.63!U:y^ 


98 


1 £+2 8 
1-6-0:5 


15.0 


< 14.40 






11.2 




r c 
JJ 


o -7+0.7 
2./-1.3 


1 n 


<14.10 






2.3 


< 14.36 






14.0 


< 14.40 






12.6 


14.60!°:^ 


425 


linear 


10.1 


18.94!0?6 


65 


3.7!?i 


15.9 


<14.10 






15.3 


< 14.26 






9.3 


< 14.20 






6.1 


18.2l!<]i^ 


41 


6.0!|t 


10.3 


19.87!;!;!!? 


91 


4 3+2.4 


27.8 


19.71!J!;S!^ 


74 


4.8!i? 


35.2 


15.55!0il 


163 


8.i!^i 


31.9 


<14.36 






6.3 




53 


2.3!i:l 


23.9 



''Color excess is total of Galactic plus LMC. Corrections are applied to correlations 
in the text. 

''LSR velocity assuming that the Galactic H2 Unes he at V/^s = km s~^ 

'^The tabulated N(H I) are corrected for the systematic effects discussed in § 2.4. 
The values here are N(H I) = (0.75 ± 0.25) N(21 cm). 
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Table 4 

FUSE H2 Small Magellanic Cloud Summary 



Name 


E(B-V)" 


vlsr'' 
(kms-1) 


N(H2) 
(cm-2) 


Ibl 

(K) 


feopp 

(kms-1) 


N(H 1)'= 
(10^" cm-2) 


AV 14 


0.135 


124±11 


17. 


m+l 35 


53 


14.3^-^ 


80 


.3 


AV 15 


0.10 


130±7 


16. 


Qr\+1 14 


47 




58 


.2 


AV26 


0.13 


122 ±9 


20. 


05 


53 




81 


.8 


AV47 


0.075 


117±8 


17. 


^^+0 34 


363 


6 1+^'' 


48, 


.0 


AV69 


0.11 


112±8 


18. 


-73+0.16 

—0.15 


90 


7 7+2.1 
' -^-2.2 


80 


.3 


AV75 


0.155 


108 ±11 


18. 


CI 40.21 
■-'^-0.62 


47 


9.6!|? 


80 


.3 


AV83 


0.12 


107 ±7 


15. 


^"-0.10 


125 


14.5!io 


27, 


.5 


AV 95 


06 

\J,\J\J 


if)7_(_7 

J- V7 / —1— / 


19. 


40+008 
^"-0.07 


78 

/ 


1 3+'-^ 

1U.J_2Q 


44 


.6 


PN2Q1 


0.105 


152 ±10 


19. 


40+0.08 
^*^-0.09 


78 




32 


.3 


NGC346-3 


0.095 


130±9 


15. 


"^-0.10 


61 


17.8^^0 


78, 


.0 


NGC346-4 


0.095 


129 ±11 


15. 


07+0.57 
■^'-0.06 


90 


19.0!3-o 


78 


.0 


NGC346-6 
HD 5980 


0.095 

0.07 


131±6 

124±10 


15. 

15. 


74+0.42 
'^-Q.IO 

'""-0.13 


105 

101 


1 7 3+5-9 


78 

26 


.0 
.1 


AV232 


0.11 


119±10 


15. 


2<;+1.38 
^"-0.19 


402 


7 0+00 
'■■-'-4.0 


27, 


.3 


Sk 82 


0.06 


139±7 


15. 


oQ+0.36 
'^y-0.18 


280 




21, 


.0 


AV238 


0.095 


151 ±13 


15. 


Qc+2.53 
^-^-0.11 


125 


20.7!5^5 


29, 


.4 


AV242 


0.085 


154 ±10 


17. 


91+1.00 
■^^-1.25 


94 


4 4+4.5 


29, 


.4 


AV264 


0.065 


170 ±22 


14. 


74+0.40 
'^-0.36 


60 


linear 


22, 


.7 


AV321 


0.055 




< 


14.56 






43 


.4 


AV327 


0.08 


110±10 


14. 


7Q+0.37 

■'='-0.26 


104 


linear 


45 


.0 


Sk 108 


0.07 




< 


14.56 






47, 


.1 


AV372 


0.115 


134±11 


16. 


cc+0.60 
,JJ_Q24 


108 


14.9!3;6 


19 


.5 


AV378 


0.085 


149 ±8 


14. 


00+0.49 
■^^-0.31 


42 


Unear 


35 


.9 


AV388 


0.105 


122 ±5 


19. 


40+0.23 
^"-0.20 


63 


3.3+^-f 


17 


.3 


Sk 159 


0.09 


141 ±11 


18. 


94+0.14 
^^-0.29 


84 




22, 


.5 


Skl88 


0.05 


118±11 


15. 


70+0.28 
'"-0.12 


91 


6.0!}:^ 


11 


.3 



^Color excess is total of Galactic plus SMC. Corrections are applied to correlations 
in the text. 

''LSR velocity assuming that the Galactic H2 Unes he at mlsr = km s"^ . 

"^The tabulated N(H I) are corrected for the systematic effects discussed in § 2.4. 
The values here are N(H I) = (0.75 ± 0.25) N(21 cm). 
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Table 5 







Rotational Level Populations - 


LMC Sight Lines 






Name 


N(0) 

(cm ^) 


N(l) N(2) N(3) 
(cm"^) (cm"^) (cm"^) 


N(4) N(5) 
(cm"^) (cm"^) 


N(6) 
(cm ^) 


^Dopp 

(kms-1) 



Sk -67 05 


19.30!«»5 


18-96-^on6 


15.66!0-32 
14.18!0-i 


15.29^^2 


14 62+0-08 








Sk -67 14 


14.18![;« 


14.86!«-j| 










linear 


HD 32109 


18.40+|}j3 


18.33+!}^^ 


15.51+,V^22 


1 5.55 A -1^ 

— U.Zd 


14.93+i-fZ 


14.58+0-68 

— u.u/ 




5.8+1-8 

—J.O 


Sk -65 21 


17.81^13 


17.94!^-^^ 


16.81!o;4o 


16.70!0-« 


14 oa+1.62 






2 3+2-1 


Sk -65 22 


14.20!«:|] 


14.52tr, 


lt.ZJ_Q J2 


14.28!o:28 








linear 


HD 33133 


17.48+''i^ 

— 1 .Ud 


16.00+H^ 

— U.iSZ 


14.74^0-76 

— U.2o 


14.50!o:50 

— U.li) 








3.5+2-5 

—1.4 


Sk -69 59 


<13.6 


<13.7 














Sk -70 60 


<13.9 


<14.1 














Sk -70 69 


<14.0 


<14.2 














Sk -68 41 


<13.8 


<13.9 














Sk -68 52 


19.28t«;^ 


19 01 -^" "^ 


15 13+' 


15 39+'-™ 


147Q+I.9O 
'^-0.21 






4 Q+2.9 


Sk -67 69 


14.82!i-5^ 


15.13^]! 


14.44!0-34 


14.74!0-85 


14.04;0;08 






6.9!^:^ 


Sk -69 104 


<13.9 


<14.6 














Sk -67 76 


14.24^0-40 


14.58+0-31 

— U.Z / 


14.33^5-^ 

—U.JO 


14.20^^-52 

— U.40 








linear 


Sk -69 124 


<13.7 


<13.9 














Sk-67 101 


<13.7 


<13.9 














Sk -67 104 


<13.6 


<13.7 














Sk -68 80 


<13.6 


<13.7 














BI 170 


<13.7 


<13.9 














Sk-67 111 


<14.3 


<14.5 














BI 173 


15.03+0-40 


ic 3C+0.53 


14.58!0-i6 


14 70+0.15 

It. /U_Q Qg 


14.20!0;06 






5 9+2-2 


Sk -70 91 


14 05+038 

lt.UJ_Q3(, 


1 4 5^+0.22 


14.28^:00 










linear 


Sk -66 100 


<13.9 


<14.0 














HD 269582 


<13.8 


<13.9 














HD 37026 




15 15+1'''' 


14 07+1.83 

^^■^'-0.09 


14 27+1 -44 








5 4+3.0 


Sk -7145 


18.06!0:[|^ 


18.26!!;;<» 


17.89!0-}7 


17.63!?:!^ 


i6.6i!?:i 


1 C QC+0.11 

1J.OJ_0.09 




1 6+2-8 
^•"-0.5 


HD 269687 


<13.9 


<14.0 














Sk -67 166 


1 5 56+i '5 


15.18!^;2o 


14.31^^^ 


14.20^:34 








9 7+0.7 
'^•'-1.3 


Sk -67 169 


<13.6 


<13.7 














Sk-67 167 


<13.8 


<14.0 














Sk -67 191 


<13.9 


<14.0 














BI 208 


13.75!«:^^ 


14.53!!J:i 












linear 


HD 37680 


18.7lt?3 


18.52^:°^ 


16.87!0-^2 


17.11!?J^ 


14.88!0-4i 


14 C!C!+0.15 




3.7!?i 


HD 269810 


<13.6 


<13.7 














BI 229 


<13.7 


<13.9 














Sk -66 169 


<13.7 


<13.8 














Sk -66 172 


18. 15^^:27 


17 •20+0.84 
^ '-^^-1.05 


1 c 4C+I.94 


15.76!^:?? 


14 50+104 






6.0^:1 


Sk -68 135 


iQ 47+0.12 


19 (:i +0.03 
^^■"^-0.03 


18.42!0:i 


18.05!":'' 


17 17+0.49 

i /-i /_1.19 


16.91+0-]^ 


14.73+0-^4 


4 3+2-4 
^■-'-2.1 


Sk -69 246 


1 Q 47+0.02 


19 18+0.03 


17 75+0.05 


17.83!0-08 


15.88!0;3^ 


15 36+0-1'* 


l4.88!o:io 


4.8!^:^ 


HD 269927 


14 71! +0.09 


1 5 23+0-30 


14 55+0.15 


14.87t38 


14.37!0:0| 






8.i!i:^ 


BI 272 


<13.8 


<14.0 














Sk -70 115 




1 Q ^7+0.07 
^^••"-0.07 


17 cc+0.34 
1 / .JJ_Q34 


17 1 9+0.14 
^ '-^^-0.81 


15.40!0;I^ 


14.80!0-« 




o + 1.4 
^••'-0.8 



The 4 a upper limits are derived assuming statistical uncertainties and a 30 km s resolution element. 



Table 6 

Rotational Level Populations'- SMC Sight Lines 



Name 


N(0) 


N(l) 


N(2) 


N(3) 


N(4) 


N(5) 


Z?Dopp 




(cm-2) 


(cm-2) 


(cm-2) 


(cm-^) 


(cm-2) 


(cm-2) 


(kms"i) 


AV 14 


16.88^^93 


16 44+1.59 


15 04+"-2'^ 


15 21+"-33 


14.65!^:^^ 




14.3!!:^ 


AV 15 


16.79!|it 


i6.i6!J:Ii 


14.93!«:g^ 


i5.i6!«:f« 


14.50!«:gf 




9.i!i:g 


AV26 


20.49!«:«^ 


20.06!«:«? 


17.99!?:^ 


i7.9i!?:«^ 


i6.74!?:^i 




i.8!?| 


AV47 


16.94!?ii 


17.69!°:2^ 




i5.30!!!:4i 


1 4 on+0.06 




6 1 +3 " 


AV69 


18.36^n;B 


18.49^«;1^^ 


15.89!>;»f, 


15.82!i:21 


14.55!!!:!!? 




7.2!|i 


AV75 


18.42!«:i 


i7.80!J!:i 


15.56!«;i 


15.72!g:« 


14.87!^:^^ 


14 75+0.10 


9.6!iY 


AV83 


14.25!"!! 


14 1+0.06 


14.49!°;!!? 


14.74!<!-!4 






linear 


AV95 


19.09tg^ 




15.30!°:?! 


i5.20!g:i: 


14.47!g:g2 




10.3!|^ 


AV207 


i9.09!«:}g 


19.09t«6 


17.47!«:i^ 


i7.i6!0:^6 


14.91^:^^ 




2.7!t:2 


NGC 346-3 


15.33!J:?« 


15.08:!fl;og 


14.66!«:i^ 


14 61+0 -15 






17.8!^7 


NGC 346-4 


14.87!«:4^ 


1 'I 00+" '''* 

iO.UU_Q07 


1 4 49+0.24 


i4.53!!!:i! 






19.0!n.3 


NGC 346-6 




1 c 44+0.52 


14 77+0.10 


14.82!!!:,',^ 






1 7 3+^-^ 


HD 5980 


15 14+«38 


iJ.JO_g [4 


i4.70!«:!S 


1 4 cQ+0.04 








AV232 


14.26!°:?^ 


10.UJ_o29 


14.38!Ji:>? 


14.51!!!:?,^ 






7 q+oo 
' .-"-4.0 


Sk 82 


14.95^16 


1 c /;4+0.48 


15.01!°:!! 


15.09!g:l| 


1 4 34+0.04 




9.4!i:^ 


AV238 


1 5 37+2-56 
IJ.J '-0.09 


1 5 73+2-60 


i4.8i!g:S 


14 72+0-49 






20.7!g.5 


AV242 


16.80!?:^« 


16.97!?i4 


15.48!j:i« 


15.27!j:f^ 






4.4!^:I 


AV264 


14 'S'S+0.42 


14.28![|;^^ 










linear 


AV321 


<14.0 


<14.2 












AV 327 


14 21+0.18 


14 4-; +0.51 


1493+0.22 








linear 


Sk 108 


<14.0 


<14.2 












AV372 


16 01+*-5^ 

1D.U1_Q jg 


16.28!°:® 


1 S 4'S+0.22 

lJ.tJ_Q J3 


15.40!^:?^ 


14 C(:+0.05 


1 4 60+" "^ 


14.9!!:^ 


AV378 


14.22!«i| 


14.48t?? 


i4.i5!«:|^ 


14.34!0:3^ 






linear 


AV388 


19.1ml 


18.96!«:i« 


i7.63:S:if 


17.24!?i? 


15.79!l:51 




^ ^+6.6 
J.J_2g 


Sk 159 


18.60!°:'^ 


18.67!«i^ 


16 60+" **" 


16 13+1'2 






4.9!?:i 


Skl88 


15.23!°^^ 


15.37!«io 


14 7(:+0.10 


14 C7+0.08 






6.0!!:^ 



The 4 a upper limits are dervied assuming statistical uncertainties and a 30 km s ' resolution element. 
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Fig. 3 . — Examples of the FUSE sample. The three stars shown are, from left to right: Sk -67 166 (also seen in Figure 4) and Sk -69 246 in the LMC and Sk 188 in the SMC. The J=0, 
1, and 2 lines of the Lyman 4-0 band are shown, plotted in velocity space. They are chosen to show the velocity separations between Galactic and Magellanic H2 and to bracket the range 
of detected H2 column densities. The left line in each pair is the Milky Way component, the right line is the Magellanic Cloud component. The signal-to-noise ratios and resolution seen 
here are typical of our sample. 
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Fig. 5. — Comparison of the gas-to-dust ratio, N(H I) / E'(B-V), where E'{B-V) for the LMC (left) and SMC (right) subtracts off the Galactic contributions. Detections of H2 are 
marked with + and non-detections are marked with X . The average ratio for the Galactic disk from Shull & van Steenberg (1985) is marked with a dotted Une. The average LMC value 
(Koornneef 1982) and the average SMC value (Fitzpatrick 1985) are marked with dashed lines. This concordance indicates that these sight lines may sample most of the diffuse neutral gas 
detected by the 21 cm emission beams, making the observed reduced molecular abundances robust. 
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Fig. 6. — Molecular fraction, /h2 = 2N(H2) / [N(H I) + 2N(H2)], versus corrected color excess. E'(B-V) was determined from FUV extinction measurements and is corrected for 
a mean Galactic E(B-V) = 0.075 for tiie LMC and E(B-V) = 0.037 for tiie SMC. For stars witii E'(B-V) < 0.075 or 0.037, we have set the corrected value to E(B-V) = 0.00. See tiie 
discussion in § 3. 
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Fig. 7. — Comparison of the FUSE LMC and SMC molecular fractions with the Copernicus (S77) sample of stars in the Galactic disk (MW). The SMC and LMC points have been 
given small random shifts in foreground-conected color excess of up to E'{B-V) = ±0.003 to separate coincident points. Below E'(B-V) = 0.08, the scatter in molecular fraction is large 
in the Clouds. For E'{B-V) > 0. 10, the SMC and LMC points lie below the average Galactic value of fni ^0.1. 
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log[N(HI)+2N(H2)] cm"^ log [N(HI) + 2N{H,)]^,, cm ' 

Fig. 8. — Left panel: Molecular fraction versus total H column density for the 44 LMC and 26 SMC sight lines and the S77 sample. Right panel: Cumulative average molecular 
fraction for the LMC, SMC, and S77 samples, as a function of the maximum value of the total gas content, N(H) = N(H I) + 2N(H2), considered in the average. The LMC and SMC are 
more gas-rich than the Galaxy, but they have a lower molecular fraction. 
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Fig. 9. — The S77, LMC, and SMC samples plotted against the total H content and compared to our numerical models. The panels are labeled with the letter of the model grids they 
display (see Table 7). Panel A: These models have Galactic values of R, 1, and the full range of kinetic temperature, gas density, and cloud size. Panel B: As before, but with R=3 X 10~'^ 
cm^ s"^ 0.1 times Galactic. Panel C: Same as before, with Galactic R and / = 10- lOOx the Galactic value. Panel D: As before, but with R = 3 X 10"^^ cm"* s"' and / = 10- lOOx the 
Galactic value. The agreement with the LMC and SMC data indicates enhanced radiation and reduced H2 formation rate in the Clouds. See § 3.2 for more discussion. The error bars are 
attached randomly to half the points to account for the systematic uncertainty in N(H I). 
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Fig. 10. — Upper left: Rotational temperature Tq\ plotted against total H2 column density N(H2). We find (7bi) = 82± 21 K f or sight lines with N{H2) > lO'^^ cm"" (excluding the 
outlier AV 47 in the SMC) and (Tqi ) = 1 15 K for all sight lines. The rotational temperatures for LMC and SMC sight lines lie near the Galactic average of 77 ± 17 K from S77, indicating 
that the physical conditions in LMC and SMC gas may be similar to those in the Galaxy if 7bi traces the kinetic temperature for high column density clouds. Upper right: The column 
density ratio N(3)/N(l) for the LMC and SMC samples, compared with the Galactic sample. There is no evidence for enhanced N(3)/N(l) in the Clouds, perhaps indicating that this ratio 
is determined by collisions at densities similar to those in Galactic gas. The blue lines are fits to the Galactic points. Lower left: The column density ratio N(4)/N(2) for the LMC and SMC 
samples, compared with the Galactic sample. Lower right: The column density ratio N(5)/N(3) for the LMC and SMC samples, compared with the Galactic sample. 
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Fig . 1 1 . — The column density ratios N(4)/N(2) and N(5)/N(3) for the LMC, SMC, and Galactic samples, compared with the model grid (in open circles). The model grids plotted have 
the range of density, size, temperature, and incident radiation field discussed in § 3.2. The LMC and SMC samples show considerable rotational excitation above the Galactic observations 
and the models, hi general, increasing incident radiation and decreasing formation rate coefficient R move model points to the lower right on this diagram. 
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Fig. 12. — IRAS 100 /iin map of the LMC (Wheelock et al. I99I), showing the survey targets. A + sign marks a detection, and X symbol designates a non-detection. The contours 
are values of 10, 100, and 1000 MJy sr~' . In the LMC, there are 23 detections of H2 in 44 sight lines. Aside from the strong detections of Ht around 30 Doradus (near 5*^ 40'", -69° ), there 
is no conelation between IRAS flux and H2 column density at a given position. The ovals mark the positions of the 22 H I supershells identified by Kim et al. (1999). The patchy nature of 
the LMC interstellar medium and superbubble blowouts may explain the low H2 detection rate in the LMC. 
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Fig. 1 3. — IRAS 100 /.im map of the SMC (Wheelock et al. 1991), showing the locations of the survey targets. A + sign marks a detection of H2, and X symbol designates a 
non-detection. The contours are values of 10, 100, and 1000 MJy sr~ . In the SMC, there are 24 detections in 26 sight lines. 
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Fig. 14. — No correlation is found between IRAS 100 /j,m surface brightness, ctfir, of the LMC CO clouds (Cohen et al. 1988) and N(H2) measured with FUSE. The method of Israel 
(1997) for deriving H2 column densities from far-infrared flux is compared to the pencil-beam FUSE survey as discussed in the text. No correlation is found between IRAS 100 pm flux 
and N(H2) measured with FUSE. The total H2 mass for the LMC extrapolated from the FUSE sample lies 10 X below the mass inferred by Israel (1997). This indicates a significantly 
lower molecular content in the LMC or a substantial cold, dense component not sampled by the FUSE sample. 



